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Abstract 

A search is made for charged Higgs bosons predicted by Two-Higgs-Doublet extensions 
of the Standard Model (2HDM) using electron-positron collision data collected by the OPAL 
experiment at ^/s = 189 — 209 GeV, corresponding to an integrated luminosity of approxi- 
mately 600 pb~^. Charged Higgs bosons are assumed to be pair-produced and to decay into 
qq, TUr or AW^*. No signal is observed. Model-independent limits on the charged Higgs- 
boson production cross section are derived by combining these results with previous searches 
at lower energies. Excluded areas on the [mjji, BR(H^— ^rz^,-)] plane are presented assum- 
ing BR(H^^Ti/,-) + BR(H^^qq) = 1. Under the above assumption, motivated by general 
2HDM type II models, charged Higgs bosons are excluded up to a mass of 76.6 GeV at 95% 
confidence level, independent of the branching ratio BR(H^^r^'T-). A scan of the 2HDM type 
I model parameter space is performed and limits on the Higgs-boson masses mjji and ttt-a 
are presented for different choices of tan /?. 
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1 Introduction 

In the Standard Model (SM) [1], the electroweak symmetry is broken via the Higgs mecha- 
nism [2] generating the masses of elementary particles. This requires the introduction of a 
complex scalar Higgs-field doublet and implies the existence of a single neutral scalar parti- 
cle, the Higgs boson. While the SM accurately describes the interactions between elementary 
particles, it leaves several fundamental questions unanswered. Therefore, it is of great interest 
to study extended models. 

The minimal extension of the SM Higgs sector required, for example, by supersymmetric 
models contains two Higgs-field doublets [3] resulting in five Higgs bosons: two charged 
(H^) and three neutral. If CP-conservation is assumed, the three neutral Higgs bosons 
are CP-eigenstates: h and H are CP-even and A is CP-odd. Two-Higgs-Doublet Models 
(2HDMs) are classified according to the Higgs-fermion coupling structure. In type I models 
(2HDM(I)) [4], all quarks and leptons couple to the same Higgs doublet, while in type II 
models (2HDM(II)) [5], down-type fermions couple to the first Higgs doublet, and up-type 
fermions to the second. 

Charged Higgs bosons are expected to be pair-produced in the process e+e~— >^H+H~ at 
LEP, the reaction c~''c~— >^H^W^ having a much lower cross section [6]. In 2HDMs, the tree- 
level cross section [7] for pair production is completely determined by the charged Higgs-boson 
mass and known SM parameters. 

The H='= branching ratios are model-dependent. In most of the 2HDM(II) parameter space, 
charged Higgs bosons decay into the heaviest kinematically allowed fermions, namely rvr and 
quark pairs-*^. The situation changes in 2HDM(I), where the decay H='=— s-AW*"* can become 

^Throughout this paper charge conjugation is imphed. For simpUcity, the notation tUt stands for r'^i'r 
and T^i>T and qq for a quark and anti-quark of any flavor. 
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dominant if the ratio of the vacuum expectation values of the two Higgs-field doublets tan /3^1 
and the A boson is sufRcicntly light [8]. 

In this paper we search for charged Higgs bosons decaying into qq, tVt and AW^* using the 
data collected by the OPAL Collaboration in 1998—2000. The results are interpreted within 
general 2HDM(II) assuming BR(H='=— >tz/^) + BR(H='=— >-qq) = 1 for the branching ratios and 
in 2HDM(I) taking into account decays of charged Higgs bosons via AW^*, as well. Our 
result is not confined to qq={cs, cs} although that is the dominant hadronic decay channel 
in most of the parameter space. 

The previously published OPAL lower limit on the charged Higgs-boson mass, under the 
assumption of BR(H='=^rz/T-) + BR(H^^qq) = 1, is mjj± > 59.5 GeV at 95% confidence level 
(CL) using data collected at < 183 GeV [9, 10]. Lower bounds of 74.4 - 79.3 GeV have 
been reported by the other LEP collaborations [11-13] based on the full LEP2 data set. The 
DELPHI Collaboration also performed a search for H^— >AW^* decay and constrained the 
charged Higgs-boson mass in 2HDM(I) [12] to be mH±>76.7 GeV at 95%CL. 

2 Experimental considerations 

The OPAL detector is described in [14]. The events are reconstructed from charged-particle 
tracks and energy deposits {clusters) in the electromagnetic and hadron calorimeters. The 
tracks and clusters must pass a set of quality requirements similar to those used in previous 
OPAL Higgs-boson searches [15]. In calculating the total visible energies and momenta of 
events and individual jets, corrections are applied to prevent double-counting of energy in the 
case of tracks and associated clusters [15]. 

The data analyzed in this paper were collected in 1998—2000 at center-of-mass energies of 
189 — 209 GeV as given in Table 1. Due to different requirements on the operational state of 
the OPAL subdetectors, the integrated luminosity of about 600 pb~^ differs slightly among 
search channels. 



Year 


1998 


1999 


2000 


Ecm (GcV) 


186-190 


190-194 


194-198 


198-201 


201-203 


200-209 


<£cm> (GeV) 


188.6 


191.6 


195.5 


199.5 


201.9 


206.0 


(GeV) 


189 


192 


196 


200 


202 


206 


JCdt (pb"') (2t) 


183.5 


29.3 


76.4 


76.6 


45.5 


212.6 


jrdt (pb-^) (2j+r,4j) 


179.6 


29.3 


76.3 


75.9 


36.6 


217.4 


(pb-^) (8j, 6j+£,4j + r) 


175.0 


28.9 


74.8 


77.2 


36.1 


211.1 



Table 1: Data-taking year, center-of-mass energy bins, luminosity-weighted average centcr-of- 
mass energies, the energies of signal and background Monte Carlo simulations, and integrated 
luminosities of the data. The data correspond to total integrated luminosities of 623.9 pb~^ 
for the two-tau, 615.1 pb~^ for the two-jet plus tau and the four-jet channels and 603.1 ph"^ 
for the H='=— >AW='=* selections. 

In this paper the following final states are sought: 

• H+H~^r+^'T-''"~^'T {two-tau final state, 2r), 

• H+H~^qqrz/T- {two-jet plus tau final state, 2j + t), 

• H"*"H~^qqqq {four-jet final state, 4j), 

• H+H-^AW+*AW-*^bbqqbbqq {eight-jet final state, 8j), 
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• H+H-^AW+*AW"*^_bbqqbb£f£ {six-jet plus lepton final state, 6j +^), 

• H+H^ — >AW^*TZ/T-^bbqqT;/T- {four-jet plus tau final state, 4j + r). 

The signal detection efficiencies and accepted background cross sections are estimated 
using a variety of Monte Carlo samples. The HZHA generator [16] is used to simulate H+H" 
production at fixed values of the charged Higgs-boson mass in steps of 1 — 5 GeV from the 
kinematic limit down to 50 GcV for fermionic decays and 40 GeV for bosonic decays. 

The background processes arc simulated primarily by the following event generators: 
PYTHIA [17] and KK2F [18] (Z/7*^qq(7)), grc4f [19] (four-fermion processes, 4f), BH- 
WIDE [20] and TEEGG [21] (e+e-(7)), KORALZ [22] and KK2F (/x+/x-(7) and t+t~{j)), 
PHOJET [23], HERWIG [24], Vermaseren [25] (hadronic and leptonic two-photon processes). 

The generated partons, both for the signal and the SM Monte Carlo simulations, are 
hadronized using JETSET [17], with parameters described in [26]. For systematic stud- 
ies, cluster fragmentation implemented in HERWIG for the process Z/7*— >^qq(7) is used. 
The predictions of 4f processes are cross-checked using EXCALIBUR [27], KoralW [28] and 
KandY [29]. 

The obtained Monte Carlo samples are processed through a full simulation of the OPAL 
detector [30] . The event selection is described below. 

3 Search for four-fermion final states 

In most of the parameter space of 2HDM(II) and with a sufficiently heavy A boson in 
2HDM(I), the fermionic decays of the charged Higgs boson dominate and lead to four-fermion 
final states. The most important decay mode is typically H^^tUt, with the hadronic mode 
H^— >qq reaching about 40% branching ratio at maximum. 

The search for the fully leptonic final state H+H~— ^-t+z/t-t^P^ is described in [31]. The 
searches for the H+H~— ^^qqri^T- and the H+H~— >^qqqq events are optimized using Monte Carlo 
simulation of H^^cs decays. The sensitivities to other quark flavors arc similar and the 
possible differences are taken into account as systematic uncertainties. Therefore, our results 
are valid for any hadronic decay of the charged Higgs boson. 

Four-fermion final states originating from H+H" production would have very similar kine- 
matic properties to W^W~ production, which therefore constitutes an irreducible background 
to our searches, especially when mH± is close to m^i. To suppress this difficult SM back- 
ground, a mass-dependent likelihood selection (similar to the technique described in [32]) is 
introduced. For each charged Higgs-boson mass tested (mtest), a specific analysis optimized 
for a reference mass (jn^^.f) close to the hypothesized value is used. 

We have chosen a set of reference charged Higgs-boson masses at which signal samples are 
generated. Around these reference points, mass regions (labeled by mref ) are defined with the 
borders centered between the neighboring points. For each individual mass region, at each 
ccntcr-of-mass energy, we create a separate likelihood selection with the signal histograms 
built using events generated at micf. The background histograms are composed of the SM 
processes and are identical for all mass regions. 

When testing the hypothesis of a signal with mass mtestj the background and data rate 
and discriminant (i.e. the reconstructed Higgs-boson mass) distribution depend on the mass 
region to which m-tcst belongs. The signal quantities depend on the value of mtcst itself and are 
determined as follows. The signal rate and discriminant distribution are computed, with the 
likelihood selection optimized for mref, for three simulated signal samples with masses miow; 
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mrcf and mhigh- Here, miow and mhigh are the closest mass points to mref at which signal Monte 
Carlo samples are generated, with miow< nT-icfK "Tihigh. The signal rate and discriminant 
distribution for rritcst arc then calculated by linear interpolation from the quantities for miow 
and mref if "T'test< ^"rcf, or for rrircf and mhigh if "T'test> rnref- 

When building the likelihood function three event classes are considered: signal, four- 
fermion background and two-fermion background. The likelihood output gives the probability 
that a given event belongs to the signal rather than to one of the two background sources. 

3.1 The two-jet plus tau final state 

The analysis closely follows our published one at s/s = 183 GeV [10]. It proceeds in two 
steps. First, events consistent with the final state topology of an isolated tau lepton, a pair of 
hadronic jets and sizable missing energy are preselected and are then processed by a likelihood 
selection. The sensitivity of the likelihood selection is improved by building mass-dependent 
discriminant functions as explained above. 

Events are selected if their likelihood output {£) is greater than a cut value chosen to 
maximize the sensitivity of the selection at each simulated charged Higgs-boson mass. Apart 
from the neighborhood of the W+W~ peak, the optimal cut does not depend significantly on 
the simulated mass and is chosen to be jC> 0.85. Around the W"'"W~ peak, it is gradually 
reduced to 0.6 at the lowest. 

At the end of the selection, 331 events are selected in the data sample with 316.9 ± 

3.2 (stat.) ±38.4 (syst.) events expected from SM processes for a test mass of mH±=75 GeV. 
The sources of systematic uncertainties are discussed below. Four-fermion processes account 
for more than 99% of the SM background and result in a large peak in the reconstructed mass 
centered at the mass (with a second peak at the Z mass for test masses of mjji > 85 GeV). 
The signal detection efficiencies for the various LEP energies are between 25% and 53% for 
any charged Higgs-boson mass. 

The likelihood output and reconstructed di-jet mass distributions for simulated Higgs- 
boson masses of 60 GeV and 75 GeV are presented in Figures l(a-d). The reconstructed 
Higgs-boson mass resolution is 2.0 — 2.5 GeV [10]. Figure 2(a) gives the mass dependence of 
the expected number of background and signal events and compares them to the observed 
number of events at each test mass. 

The systematic uncertainties are estimated for several choices of the charged Higgs-boson 
mass from 50 GeV to 90 GeV at ccnter-of-mass energies of ■y/i=189 GeV, 200 GeV and 
206 GeV to cover the full LEP2 range. The following sources of uncertainties are considered: 
limited number of generated Monte Carlo events, statistical and systematic uncertainty on 
the luminosity measurement, modeling of kinematic variables in the pre- and likelihood se- 
lections, tau lepton identification, dependence of the signal detection efficiency on final-state 
quark flavor, signal selection efficiency interpolation between generated Monte Carlo points, 
background hadronization model, and four-fermion background model. The contributions 
from the different sources are summarized in Table 2. 

In the limit calculation, the efficiency and background estimates of the 2] + t channel 
are reduced by 0.8—1.7% (depending on the center-of-mass energy) in order to account for 
accidental vetoes due to accelerator-related backgrounds in the forward detectors. 
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OPAL: H*H"^qqTV, 189-209 GeV 



OPAL: H*H"^qqTV, 189-209 GeV 



± = 60 GeV 
OPAL data 
I I signal C60 GeV 




O 0.2 0.4 0.6 0.8 1 
Likelihood 
OPAL: H*H"^qqTV, 189-209 GeV 




60 80 100 

Reconstructed mass (GeV) 
OPAL: H*H"^qqTV, 189-209 GeV 




O 0.2 0.4 0.6 0.8 1 
Likeliliood 

OPAL: H*H"^qqqq, 189-209 GeV 

. 31 



60 80 100 

Reconstructed mass (GeV) 



OPAL: H H'^qqqq, 189-209 GeV 



„± = 60GevV£l J 

OPAL data ^ ' E 
I I signal (60 GeV) 

I 1 4f 

H 2f 




0.2 0.4 0.6 0.8 1 
Likeliliood 
OPAL: H*H"^qqqq, 189-209 GeV 



80 100 
Reconstructed mass (GeV) 
OPAL: H^H'^qqqq, 189-209 GeV 




0.2 0.4 0.6 0.8 1 
Likelihood 



60 80 100 

Reconstructed mass (GeV) 



Figure 1: Likelihood output and reconstructed di-jet mass distributions for the (a-d) 2 j + r 
and (e-h) 4j channels. The distributions are summed up for all center-of-mass energies and 
correspond to 60 GeV and 75 GeV simulated charged Higgs-boson masses. All Monte Carlo 
distributions are normalized to the integrated luminosity of the data. When plotting the 
likelihood output, the signal expectation is scaled up by a factor of 10 for better visibility. A 
hadronic branching ratio of 0.5 is assumed for the 2 j + r signal, and 1.0 for the 4j signal. The 
reconstructed mass distributions are shown after the likelihood selection. 
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OPAL: H*H qqxv, 189-209 GeV OPAL: H*H qqqq, 189-209 GeV 




50 60 70 80 90 100 50 60 70 80 90 



m„± (GeV) m„± (GeV) 

Figure 2: The number of observed data, expected background and signal events for the (a) 
2 j + r and (b) 4j channels. The numbers are summed up for all center-of-mass energies and 
shown as a function of the reference charged Higgs-boson mass. A hadronic branching ratio 
of 0.5 is assumed for the 2j + r signal, and 1.0 for the 4j signal. Each bin corresponds to a 
different likelihood selection optimized for the mass at which the dot is centered. Since the 
same background simulations are used to form the reference histograms and the same data 
enter the selection, the neighboring points are strongly correlated. 



Source 


2j+r 


4j 


signal 


background 


signal 


background 


MC statistics 


3.1-4.6 


1.4-4.3 


1.6-2.4 


0.9-1.9 


luminosity 


0.3 


0.3 


0.3 


0.3 


preselection 


1.5-4.7 


1.8-7.6 


0.3-1.1 


0.5-2.2 


likelihood selection 


0.9-6.5 


5.8-22.7 


0.7-2.4 


2.1-7.5 


tau identification 


3.0 


3.0 


N.A. 


N.A. 


quark flavor 


2.7-3.8 


N.A. 


1.2-6.4 


N.A. 


interpolation 


0.2-0.4 


N.A. 


0.7-3.7 


N.A. 


hadronization model 


N.A. 


1.0-2.7 


N.A. 


0.7-4.1 


4f background model 


N.A. 


0.3-3.3 


N.A. 


1.7-3.7 



Table 2: Relative systematic uncertainties on the expected background and signal rates for 
the 2j + r and 4j Gnal states. The numbers are given in % and depend on the centre-of-mass 
energy and the reference charged Higgs-boson mass. N.A. stands for not applicable. 
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3.2 The four-jet final state 

The event selection follows our published analysis at i/i=183 GeV [10]: first, well-separated 
four-jet events with large visible energy are preselected; then a set of variables is combined 
using a likelihood technique. To improve the discriminating power of the likelihood selec- 
tion, a new reference variable is introduced: the logarithm of the matrix element probability 
for W+W~ production averaged over all possible jet-parton assignments computed by EX- 
CALIBUR [27]. Moreover, we introduce mass-dependent likelihood functions as explained 
above. As the optimal cut value on the likelihood output is not that sensitive to the charged 
Higgs-boson mass in this search channel, we use the condition L > 0.45 at all center-of-mass 
energies and for all test masses. 

There is a good agreement between the observed data and the SM Monte Carlo expec- 
tations at all stages of the selection. After all cuts, 1100 events are selected in the data, 
while 1117.8 ± 5.9 (stat.) ±74.4 (syst.) events are expected from SM processes for a test 
mass of mH±=75 GeV. The four-fermion processes account for about 90% of the expected 
background and result in a large peak centered at the mass and a smaller one at the Z 
boson mass. The signal detection efficiencies are between 41% and 59% for any test mass and 
center-of-mass energy. 

Typical likelihood output and reconstructed di-jet mass distributions of the selected events 
together with the SM background expectation and signal shapes for simulated charged Higgs- 
boson masses of 60 GeV and 75 GeV are plotted in Figures l(e-h). The Higgs-boson mass 
can be reconstructed with a resolution of 1 — 1.5 GeV [10]. Figure 2(b) shows the mass 
dependence of the expected number of background and signal events and compares them to 
the observed number of events at each test mass. Systematic uncertainties are estimated in 
the same manner as for the 2j +t search and are given in Table 2. 

4 Search for AW+*AW~* events 

In a large part of the 2HDM(I) parameter space, the branching ratio of H^ — ^ AW^* dom- 
inates. The possible decay modes of the A boson and the W^* lead to many possible 
H+H~^AW"'"*AW~* event topologies. Above mA~12 GeV, the A boson decays predomi- 
nantly into a bb pair, and thus its detection is based on b-flavor identification. Two pos- 
sibilities, covering 90% of the decays of two W='=*, arc considered: quark pairs from both 
W^* bosons or a quark pair from one and a Icptonic final state from the other. The event 
topologies are therefore "eight jets" or "six jets and a lepton with missing energy", with four 
jets containing b-fiavor in both cases. 

The background comes from several Standard Model processes. ZZ and W+W" produc- 
tion can result in multi-jet events. While ZZ events can contain true b-flavored jets, W^W~ 
events are selected as candidates when c-flavored jets fake b-jets. Radiative QCD corrections 
to e"'"e~-^qq also give a significant contribution to the expected background. 

Due to the complexity of the eight-parton final state, it is more efficient to use general event 
properties and variables designed specifically to discriminate against the main background 
than a full reconstruction of the event. As a consequence, no attempt is made to reconstruct 
the charged Higgs-boson mass. 

The analysis proceeds in two steps. First a preselection is applied to select b-tagged multi- 
jet events compatible with the signal hypothesis. Then a likelihood selection (with three event 
classes: signal, four-fermion background and two-fermion background) is applied. 
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The preselection of multi-jet events uses the same variables as the search for the hadronic 
final state in [10] with optimized cut positions. However, it introduces a very powerful new 
criterion, especially against the W+W" background, on a combined b-tagging variable (^Sevt) 
requiring the consistency of the event with the presence of b-quark jets. 

The neural network method used for b-tagging in the OPAL SM Higgs-boson search [15] 
is used to calculate on a jet-by-jet basis the discriminating variables f^^^ and /^^g/^- These 
are constructed for each jet i as the ratios of probabilities for the jet to be c- or uds-like 
versus the probability to be b-like. The inputs to the neural network include information 
about the presence of secondary vertices in a jet, the jet shape, and the presence of leptons 
with large transverse momentum. The Monte Carlo description of the neural network output 
was checked with LEPl data with a jet energy of about 46 GeV. The main background in 
this search at LEP2 comes from four-fermion processes, in which the mean jet energy is about 
50 GeV, very close to the LEPl jet energy; therefore, an adequate modeling of the data is 
expected with an event reconstruction assuming four jets. 

The AW+*AW~* signal topology depends on the Higgs-boson masses. At mA~12 GeV 
or mA«mH±, the available energy in the A or W^* system is too low to form two clean, 
coUimated jets. At high mH± , the boost of the A and W='=* bosons is small in the laboratory 
frame and the original eight partons cannot be identified. At low mjji, the A and W^* 
bosons might have a boost, but it is still not possible to resolve correctly the two partons 
from their decay. Prom these considerations, one can conclude that it is not useful to require 
eight (or even six) jets in the event, as these jets will not correspond to the original partons. 
Consequently, to get the best possible modeling of the data, four jets are reconstructed with 
the Durham jet-finding algorithm [33] before the b-tagger is run. 

The flavor-discriminating variables are combined for the four reconstructed jets by 

The index i runs over the reconstructed jets (f = 1, ...4) and the parameters a and (3 are 
numerical coefficients whose optimal values depend on the flavor composition of the signal 
and background final states. However, since the expected sensitivity of the search is only 
slightly dependent on the values of a and (3, they are fixed at a = 0.1 and /3 = 0.7. Events 
are retained if Bevt > 0.4. 

The preselections of the two event topologies (8j and 6j + i) arc very similar. However, in 
the 6i+£ channel, no kinematic fit is made to the W+W^-^qqqq hypothesis and, therefore, 
no cuts are made on the fit probabilities. No lepton identification is applied; instead the 
search is based on indirect detection of the associated neutrino by measuring the missing 
energy. 

After the preselection the observed data show an excess over the predicted Monte Carlo 
background. This can partly be explained by the apparent difference between the gluon 
splitting rate into cc and bb pairs in the data and in the background Monte Carlo simulation. 
The measured rates are g^c = 3.2 ± 0.21 ± 0.38% [34] and ^^B = 0-307 ± 0.053 ± 0.097% [35] 
from the LEPl OPAL data. The gluon splitting rates in the Monte Carlo are extracted from 
e+e-^ZZ^£+^-qq events and are found to be g^'^ = 1.33 ± 0.06% and g^ = 0.116 ± 
0.0167%, averaged over all center-of-mass energies. This mismodeling can be compensated 
by reweighting the SM Monte Carlo events with gluon splitting to heavy quarks by universal 
reweighting factors [36] and at the same time deweighting the non-split events to keep the 
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total numbers of W+W", ZZ and two-fermion background events fixed at generator level. 
The reweighting factor is 2.41 for g^cc and 2.65 for g^bb. This correction results in a 
background enhancement factor of 1.08 to 1.1 after the preselection, depending on the search 
channel and the center-of-mass energy, but it does not affect the shape of the background 
distributions. The numbers of preselected events after the reweighting are given in Table 3. 



LEP energy 


preselection 


preselection 


exclusive 


exclusive 


overlap 


(year) 


8j 


Gj + (. 


8j 


Gj + /' 




189 GeV data 


238 


358 


3 


24 


5 


(1998) background 


231.2±2.9 


342.2±3.6 


2.1±0.3 


24.4±1.0 


6.3±0.5 


192 - 202 GeV data 


297 


310 


16 


16 


17 


(1999) background 


270.4±2.9 


285.0±3.0 


13.3±0.7 


10.4±0.6 


13.4±0.7 


200 - 209 GeV data 


265 


281 


9 


8 


15 


(2000) background 


252.5±3.7 


270.5±5.0 


13.0±0.9 


9.3±0.8 


12.9±0.9 



Table 3: Observed data and expected SM background events for each year in the AW"'"*AW~* 
searches. The 8j and 6 j + £ event samples after the preselection step (2nd and 3rd columns) 
are highly overlapping. After the likelihood selection, the overlapping events are removed 
from the 8j and 6j + ^ samples and form a separate search channel (last three columns). The 
uncertainty on the background prediction due to the limited number of simulated events is 
given. The Monte Carlo reweighting to the measured gluon splitting rates is included. 

As a final selection, likelihood functions are built to identify signal events. The reference 
distributions depend on the LEP energy, but they are constructed to be independent of the 
considered (mH± ) "^a) combination. To this end, we form the signal reference distributions 
by averaging all simulated H+H" samples in the {m}i±,mA) mass range of interest. 

Since the selections at ^/s = 192 — 209 GeV are aimed at charged Higgs-boson masses 
around the expected sensitivity reach of about 80—90 GeV, all masses up to the kinematic 
limit are included. On the other hand, at v^=189 GeV only charged Higgs-boson masses 
up to 50 GeV are included since the selections at this energy are optimized to reach down 
to as low as a charged Higgs-boson mass of 40 GeV where the LEPl exclusion limit lies. 
The input variables for the 8j final state are: the Durham jet-resolution parameters^ logio 2/34 
and logioy56) the oblateness [37] event shape variable, the opening angle of the widest jet 
defined by the size of the cone containing 68% of the total jet energy, the charge-signed 
cosine of the production angle in the W^W"— ^qqqq hypothesis, and the b-tagging variable 
Bcvt- At -^5=189 GeV, log^gyss, logioy45, logioye?) and the maximum jet energy are also 
used. Moreover, the sphericity [38] event shape variable has more discriminating power and 
thus replaces oblateness. Although the yij variables are somewhat correlated, they contain 
additional information: their differences reflect the kinematics of the initial partons. 

The input variables for the 6] + i selection are: log 2/34, log^Qj/sg, the oblateness, the 
missing energy of the event, and i3evt- At y^=189 GeV, logj^gysa, the maximum jet energy 
and the sphericity are also included. 

Events are selected if they pass a lower cut on the likelihood output. This final cut 
docs not rcmow (especially in the year 1999 data) all the excess events observed after the 

"^Throughout this paper yij denotes the parameter of the Durham jet finder at which the event classification 
changes from i-jet to j-jet, where j = i + 1. 
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preselection. The distributions of the critical selection variables, 2/34, y^Q and i3evt, are plotted 
on Figure 3 both in background-enriched data samples and after the preselection. To prepare 
these background-enriched data samples, the preselection cuts on 1/34 and B^yt dropped, 
except for the study of the yse variable where we keep the cut on 1/34 in order to select multi-jet 
events. The resulting samples are completely dominated by background processes. We find 
systematic differences between the data samples and the equivalent background Monte Carlo 
simulations at both stages. The observed excess cannot be attributed to a Higgs-boson signal 
whose contribution is expected to be at most 0.5% in the background-enriched data samples. 
The interpretation of the excess in terms of a systematic uncertainty on the background 
prediction will be discussed later. 

The non-signal nature of the excess is also demonstrated by Figure 4, which shows the 
distributions of the likelihood output. The positions of the likelihood cut are indicated by 
vertical lines. 

Because the 8j and the 6] + i selections have several discriminating variables in common, 
a major overlap between these selections exists. To assure that every event is counted only 
once, the two samples are redistributed into three: (i) events exclusively classified as 8j 
candidates, (m) events exclusively classified as 6j -I- ^ candidates and (in) events accepted 
by both selections. If an event falls into class (Hi), the larger likelihood output of the two 
selections is kept for further processing. The final results using the above classification are 
quoted in Table 3. 

This modified channel definition not only removes the overlap but also increases the ef- 
ficiency for detecting signal events by considering the cross-channel efficiencies (e.g. the 
efficiency to select H"'"H^^bbqqbbqq signal by the exclusive 6] + i selection can be as high 
as 18%, though it is typically only a few %). The efficiencies are determined independently 
for all simulated (mH±,wiA) combinations and interpolated to arbitrary (mH±,mA) by two- 
dimensional spline interpolation. The behavior of the selection efficiencies depends strongly 
on the targeted charged Higgs-boson mass range and also varies with the mass difference 
Am = mH± — ruA- In most cases the overlap channel has the highest efficiency. At y^=189 
GeV and mH±=45 GeV, it reaches 32% for the bbqqbb^f^ and 44% for the H+H~— ^^bbqqbbqq 
signal close to the mjj±=mA diagonal. At ■y/s=206 GeV and mH±=90 GeV, the overlap effi- 
ciency can be as high as 62% for the bbqqbb^f^ and 71% for the H+H~-^bbqqbbqq signal. 
The exclusive 6j -I- ^ selection has efficiencies typically below 20—30%, while the exclusive 8j 
selection below 10—15%. Table 4 gives the selection efficiencies at selected {m^±,mA) points. 

The composition of the background depends on the targeted Higgs-boson mass region. 
In the low-mass selection (■y/i=189 GeV) that is optimized for mjj±=40— 50 GeV, the Higgs 
bosons are boosted and therefore the final state is two-jet-like with the largest background 
contribution coming from two-fermion processes: they account for 52% in the exclusive 8j, 
80% in the exclusive + i and 76% in the overlap channel. On the other hand, in the high- 
mass analysis (-^5 = 192 — 209 GeV) the four-fermion fraction is dominant: it is 69% in the 
8j, 56% in the 6 j + ^ and 70% in the overlap channel. 

Although our studies show that the observed excess does not lie preferentially in the 
phase space region where the Higgs-boson signal is expected and that it probably originates 
from deficiencies of the Monte Carlo description of the jet-resolution parameters yij and the 
b-tagging variable B^vt, no further correction is applied to the estimated background in the 
background subtraction procedure of the statistical analysis to calculate exclusion limits. This 
leads to conservative limits on the production cross section. Therefore, no systematic uncer- 
tainty is assigned to the modeling of yu and Bevt at the preselection level. The systematic 
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OPAL: H^H^bbqqbbqq, 192-209 GeV 
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Figure 3: Most important selection variables: (a-b) logxoy34, (c-d) log^o 2/56 and (e-f) Bevt in 
the 8j channel at ^/s = 192 — 209 GeV. The distributions are shown (left) in a background- 
enriched data sample (see text for explanation) and (right) after the full preselection. To 
form the signal histograms, the Monte Carlo distributions are averaged for all simulated 
(mjj±,mA) mass combinations in the mass range of interest. The Monte Carlo reweighting 
to the measured gluon splitting rates is included. The expectations from SM processes are 
normalized to the data luminosity. The preselection cuts on 2/34 and Be^t are indicated by 
vertical lines. 13 




OPAL: ff'H^ bbqqbblv, 189 GeV OPAL: H^ff^ bbqqbblv, 192-209 GeV 




Likelihood Likelitiood 



Figure 4: Likelihood output distributions for the (a-b) 8j and (c-d) 6] + i channels at 
y/s=189 GeV and 192 — 209 GeV. To form the signal histograms, the Monte Carlo dis- 
tributions are averaged for all simulated (mH±,mA) mass combinations in the mass range of 
interest. The Monte Carlo reweighting to the measured gluon splitting rates is included. The 
expectations from SM processes are normalized to the data luminosity. The lower cuts on the 
likelihood output are indicated by vertical lines. 
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(mH±,mA) 


(GcV, GcV) 






(45,30) 


(80,50) 


(45,30) 


(90,60) 






^=189 GeV 


^=206 GcV 


bbqqbbqq 


8j 


4.6 


1.0 


9.3 


12.4 




overlap 


41.0 


2.8 


14.9 


69.6 




6j+£ 


17.0 


3.6 


4.1 


3.1 




total 


62.6 


7.4 


28.3 


85.0 


bbqqbb^t'^ 


6j+£ 


28.2 


6.0 


11.6 


7.0 




overlap 


31.8 


3.6 


14.2 


62.2 




8j 


1.8 


0.1 


2.2 


6.1 




total 


61.8 


9.7 


28.0 


75.3 


bbqqri^T- 


4j + r 


68.0 


0.0 


12.3 


11.1 



Table 4: Signal selection efRciencies in percent for the H^^AW^* final states in the different 
search channels at ^/s=189 and 206 GeV at representative (mH±,?TiA) points. 

uncertainties related to the other preselection variables, estimated from background-enriched 
data samples, are taken into account. The effect of mismodeling the shapes of the reference 
distributions in the likelihood selection is estimated for all variables, including and Bevt 
and is accounted for in the statistical analysis. 

Systematic uncertainties arise also due to the gluon splitting correction: the experimental 
uncertainty on the gluon splitting rate translates into uncertainties on the total background 
rates. Moreover, the weighted background counting introduces an uncertainty due to the 
Monte Carlo statistics of the g^cc and bb events. 

In summary, the following sources of systematic uncertainties are considered: limited 
number of simulated signal and background events, modeling of the preselection variables 
(other than y^j and Bg^t), modeling of the shapes of the reference distributions in the likelihood 
selection and the gluon splitting correction. Uncertainties below the 1% level are neglected. 
The different contributions are summarized in Table 5. 

5 Search for AW^*TUr events 

In some parts of the 2HDM(I) parameter space, both the fermionic H^^rz/,- and the bosonic 
11=*=— >AW^* decay modes contribute. To cover this transition region parallel to the (mH± , mA) 
diagonal, a search for the final state H+H~— >^AW='=*tz/t- is performed. The transition region 
is wide for small tan/3 and narrow for large tan/3; therefore, this analysis is more relevant for 
lower values of tan/3. 

Only the hadronic decays of W^* and the decay A^bb are considered. Thus the events 
contain a tau lepton, four jets (two of which are b-flavored) and missing energy. Separating 
the signal from the W^W~ background becomes difficult close to m}i±=m-\^±. 

The preselection is designed to identify hadronic events containing a tau lepton plus 
significant missing energy and transverse momentum from the undetected neutrino. In most 
cases it is not practical to reconstruct the four jets originating from the AW='=* system. Instead, 
to suppress the main background from semi-leptonic W"'"W~ events, we remove the decay 
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Source 


exclusive 8j 


exclusive 6j + ^ 


overlap 


signal 


background 


signal 


background 


signal 


background 


MC statistics 


>15/>8.5 


13.2/6.7 


>5.7/>8.4 


4.0/8.3 


>4.0/>2.8 


7.9/7.1 


preselection 


1.0 


1.0 


1.0 


2.0 


1.0 


2.0/1.5 


L selection 
















4.0/1.8 


6.0/6.2 


2.2/2.6 


6.0/8.0 


2.8/1.5 


5.5/4.9 


b-tag 


0.0/1.8 


4.7/7.0 


2.9/1.4 


4.4/7.1 


1.0/1.3 


4.0/5.1 


other 


1.8/0.7 


3.9/3.2 


1.0/1.6 


3.7/3.5 


1.2/0.7 


3.4/2.4 


gluon splitting 














g— >cc, exp. 


N.A. 


0.6/2.5 


N.A. 


1.6 


N.A. 


1.4/2.3 


g^cc, MC 


N.A. 


0.2/0.8 


N.A. 


0.6 


N.A. 


0.5/0.8 


g^bb, exp. 


N.A. 


1.4/4.2 


N.A. 


3.8/4.3 


N.A. 


5.5/5.4 


g^bb, MC 


N.A. 


0.5/1.7 


N.A. 


1.5/1.7 


N.A. 


2.2 


gluon splitting 
correction factor 


N.A. 


1.05/1.18 


N.A. 


1.13/1.15 


N.A. 


1.15/1.22 



Table 5: Relative systematic uncertainties in percent for the AW'''*AW~* searches. Where 
two values are given separated by a "/", the first belongs to the 189 GeV selection and the 
second to the 192 — 209 GeV selections. For the signal, the uncertainties due to the limited 
Monte Carlo statistics are calculated by binomial statistics for a sample size of 500 events 
and they also depend, via the selection efficiency, on the assumed Higgs-boson masses. N.A. 
stands for not applicable. The multiplicative gluon splitting correction factors, used to obtain 
the background-rate estimates as explained in the text, are given in the last line. 

products of the tau candidate and force the remaining hadronic system into two jets by the 
Durham algorithm. The requirements are then based on the preselection of Section 3.1 with 
additional preselection cuts on the effective center-of-mass energy, logj^Q j/12 and logj^g 2/23 of 
the hadronic system, and the charge-signed production angle. 

The Hkehhood selection uses seven variables: the momentum of the tau candidate, the 
cosine of the angle between the tau momentum and the nearest jet, log^^Q yi2 of the hadronic 
system, the cosine of the angle between the two hadronic jets, the charge-signed cosine of 
the W='= production angle, the invariant mass of the hadronic system, and the b-tagging 
variable Bc^t- Here, Bcvt is defined using the two jets of the hadronic system using Eq. (1) 
of Section 4, with i = 1,2 and a = f3 = I. To form the signal reference distributions, all 
simulated H'^H" samples in the (mH± , wa) mass range of interest are summed up. Since the 
search at i/s = 192—209 GeV targets intermediate charged Higgs-boson masses (60—80 GeV), 
all masses up to the kinematic limit are included. At v^=189 GeV, only charged Higgs-boson 
masses up to 50 GeV are included since the selection is optimized for low charged Higgs-boson 
masses (40-50 GeV). 

The likelihood output distributions are shown in Figure 5. There is an overall agreement 
between data and background distributions, apart from a small discrepancy at 189 GeV. 
Events are selected if their likelihood output is larger than 0.9. In total, 15 data events survive 
the selection at ^/s = 192 - 209 GeV, to be compared with 14.8 ± 0.6 (stat.) ±1.9 (syst.) 
events expected from background sources. At y^=189 GeV, where the selection is optimized 
for low Higgs-boson masses, 13 data events are selected with 6.1 it 0.5 (stat.) ±1.3 (syst.) 
events expected. The contribution of four-fermion events, predominantly from semi-leptonic 
W+W" production, amounts to 33% at ^^=189 GeV and to 90% at ^/s = 192 - 209 GeV. 
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Likelihood Likeliliood 

Figure 5: Likelihood output distribution for the 4j + r channel at (a) y/s=189 GeV and 
(b) ^/s=192 — 209 GeV. To form the signal histograms, the Monte Carlo distributions are 
averaged for all simulated {m^± , mA) mass combinations in the mass range of interest. The 
expectations from SM processes are normalized to the data luminosity. The lower cuts on the 
likelihood output are indicated by vertical line. 

At ^/s = 192 — 209 GeV, the signal selection efficiency starts at about 5% at mH±=40 GeV, 
reaches its maximum of about 40% (depending on the mass difference Am = — m.^) 
at mH±=60 GeV, then decreases to 12% at mH±=90 GeV. In the low- mass selection at 
■^5=189 GeV, the efficiency depends strongly on the mass difference: at mH±=40 GeV, it is 
22% for Am = 2 GeV and 60% for Am = 10 GeV. The selection efficiency approaches its 
maximum at mH±=50 GeV (70% for Am = 10 GeV) and then drops to zero at mjj±=80 GeV. 
Table 4 gives selection efficiencies at representative {mii±,mA) points. 

The systematic uncertainties due to the modeling of selection variables are evaluated with 
the method developed for the AW+*AW~* channels and summarized in Table 6. 



Source 


4j +r 




signal 


background 


MC statistics 


>2.7/>4.5 


8.2/7.0 


preselection: tau ID 


0.0 


2.3/5.0 


other 


0.0/1.0 


9.5/7.9 


likelihood selection: b-tag 


0.3/1.4 


2.4/3.0 


other 


3.2/2.1 


18.4/11.2 



Table 6: Systematic uncertainties in percent for the 4j + r channel. Where two values are 
given separated by a "/", the first one belongs to the 189 GeV selection and the second to the 
192 — 209 GeV selections. For the signal, the uncertainties due to the limited Monte Carlo 
statistics are calculated by binomial statistics for a sample size of 500 events and they also 
depend, via the selection efRciency, on the assumed Higgs-boson masses. 
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6 Interpretation 



None of the searches has revealed a signal-hke excess over the SM expectation. The re- 
sults presented here and those published previously [10, 31] by the OPAL Collaboration are 
combined using the method of [39] to study the compatibility of the observed events with 
"background-only" and "signal plus background" hypotheses and to derive limits on charged 
Higgs-boson production. The statistical analysis is based on weighted event counting, with 
the weights computed from physical observables, also called discriminating variables of the 
candidate events (see Table 7). Systematic uncertainties with correlations are taken into ac- 
count in the confidence level (p- value) calculations. To improve the sensitivity of the analysis, 
they are also incorporated into the weight definition [39]. 



Channel 


(GeV) 


Discriminant 


2r 


183 


simple event counting 


2r 


189-209 


likelihood output 


2j + r 


183-209 


reconstructed di-jet mass 


4j 


183-209 


reconstructed di-jet mass 


8j 


189 


simple event counting 


8j 


192-209 


likelihood output 


6j + £ 


189 


simple event counting 


6j + £ 


192-209 


likelihood output 


4j+T 


189-209 


simple event counting 



Tabic 7: Discriminating variables entering the statistical analysis for each search topology. 
Previously published results are also included. 

The results are interpreted in two different scenarios: in the traditional, super symmetry- 
favored 2HDM(II) (assuming that there are no new additional light particles other than the 
Higgs bosons) and in the 2HDM(I) where under certain conditions fermionic couplings are 
suppressed. 

First, we calculate 1 — CLf,, the confidence [39] under the background-only hypothesis, and 
then proceed to calculate limits on the charged Higgs-boson production cross section in the 
signal -I- background hypothesis. These results are used to provide exclusions in the model 
parameter space, and in particular, on the charged Higgs-boson mass. 

2HDM Type II 

First a general 2HDM(II) is considered, where BR(H^^rz/^) -|- BR(H^^qq) = 1. This model 
was thoroughly studied at LEP. It is realized in super symmetric extensions of the SM if no 
new additional light particles other than the Higgs bosons are present. As our previously 
published mass limit in such a model is mH± > 59.5 GeV [10], only charged Higgs-boson 
masses above 50 GcV are tested. Cross-section limits for lower masses can be found in [9]. In 
this model, the results of the 2r, 2] + t and 4j searches enter the statistical combination. 

The confidence 1 — CLh is plotted for each channel separately in Figure 6(a) and com- 
bined in Figure 6(b). Note that 1 — CLj, < 0.5 translates to negative values of sigma (as 
indicated by the dual y-axis scales in Figure 6(a)) and indicates an excess of events. The 
largest deviation, a 2.8a excess observed in the 2j + t channel at mH±=88 GeV, comes 
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from a known deficiency [40] of the Monte Carlo simulation of isolated tracks from the frag- 
mentation and hadronization process. When the three final states are combined assuming 
BR(H^— >tz/t-) + BR(H^— >qq) = 1, no deviation reaches the 2a level. 




Figure 6: The observed confidence levels for the background interpretation of the data, 1 — 
CLb, (a) for the three different final states as a function of the charged Higgs-boson mass, and 
(b) for the combined result in 2HDM(II) assuming BR(H=*=^TZ^r) + BR(H±^qq) = 1 on the 
[m,jj±, BR(H^^rz^T-)] plane. The significance values corresponding to the different shadings 
are shown by the bar at the right. 

The results are used to set upper bounds on the charged Higgs-boson pair production 
cross section relative to the 2HDM prediction as calculated by HZHA. The limits obtained 
are shown for each channel separately in Figures 7(a-c) and combined in Figure 7(d). The 
combined results are shown by "isolines" along which iT95(H+H~)/(T2hdm, the ratio of the 
limit on the production cross section and the 2HDM cross-section prediction, is equal to the 
number indicated next to the curves. 

Excluded areas on the [mjj± , BR(H^^r^'T-)] plane are presented for each channel sepa- 
rately in Figure 8(a) and combined in Figure 8(b). The expected mass limit from simulated 
background experiments, assuming no signal, is also shown. For the combined results, the 
90% and 99% CL contours are also given. Charged Higgs bosons are excluded up to a mass 
of 76.6 GeV at 95% CL, independent of BR(H±^rz/^). Lower mass limits for different values 
of BR(H^— >Ti/^) are presented in Table 8. 
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OPAL: H"^ H" ^ TVTV OPAL: H"^ H" qqxv 




mjj± (GeV) mjj± (GeV) 

Figure 7: Observed and expected 95% CL upper limits on the H+H^ production cross section 
times the relevant decay branching ratios relative to the theoretical prediction for the (a) 
rvT-rVr (b) qqTiy-r and (c) qqqq channels. The horizontal lines indicate the maximum possible 
branching ratios for a given channel. In (b), BKqqri/ = 2 • BR(H^^t^'t-) • BR(H^— >qq). 
(d) Upper limits on the production cross section relative to the 2HDM prediction on the 
[mH±, BR(H=''^rz^T-)] plane in 2HDM(n) assuming BR(H^^TUr) + BR(H±^qq) = 1. The 
plotted curves are isolines along which the limit is equal to the number indicated. 
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Figure 8: Observed and expected excluded areas at 95% CL on the [mH± , BR(H^— ^ri/,-)] 
plane (a) for each search channel separately and (h) combined in 2HDM(n) assuming 
BR(H=^^rz^^) + BR(H=^^qq) = 1. For the combined result, the 90% and 99% CL observed 
limits are also shown. See Table 8 for numerical values of the combined limit. 





Lower mass 


limit (GeV) 


Observed 


Expected 





76.9 


77.9 


0.5 


79.2 


78.0 


0.65 


82.0 


81.7 


1 


91.2 


89.2 


any 


76.6 (0.15) 


76.8 (0.27) 



Table 8: Observed and expected lower limits at 95% CL on the mass of the charged Higgs 
boson in 2HDM(II) assuming liI{(H^^Tiyr) + BR(H^^qq) = 1. For the results independent 
of the branching ratio (last line), the BR(H^— ^rz^,-) vaiue at which the limit is set, is given in 
parenthesis. 
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2HDM Type I 

We present here for the first time an interpretation of the OPAL charged Higgs-boson searches 
in an alternative theoretical scenario, a 2HDM(I). The novel feature of this model with respect 
to the more frequently studied 2HDM(II) is that the fermionic decays of the charged Higgs 
boson can be suppressed. If the A boson is light, the H^^AW^* decay may play a crucial 
role. 

The charged Higgs-boson sector in these models is described by three parameters: mjji, 
ruA and tan/3. To test this scenario, the charged Higgs-boson decay branching ratios are 
calculated by the program of Akeroyd et al. [8], and the model parameters are scanned in 
the range: 40 GeV < mH± < 94 GeV, 12 GeV < niA < mH±, 0.1 < tan/3 < 100. Charged 
Higgs-boson pair production is excluded below 40 GeV by the measurement of the Z boson 
width [41]. As the A boson detection is based on the identification of b-quark jets, no limits 
are derived for mA< 2m\). Below tan/3=0.1, BR(H^^AW^*) vanishes and the limit is no 
longer sensitive to ruA- 

Both the fermionic (2r, 2j + t and 4j) and the bosonic (4j +t, 6i + £ and 8j) final states 
play an important role and therefore their results are combined. There is, however, a sig- 
nificant overlap between the events selected by the H'^H"— >qqqq and H+H^^AW+*AW~* 
selections, and the events selected by the H+H~— ^qqri/',- and H+H^— ^AW='=*r^',- selections. 
Therefore, an automatic procedure is implemented to switch off the less sensitive of the over- 
lapping channels, based on the calculation of the expected limit assuming no signal. In general 
the fermionic channels are used close to the [m^±,mA] diagonal and for low tan/3, and the 
searches for H='=— >AW^* are crucial for low values of niA and high values of tan /3. 

The confidence 1 — CL^ is calculated combining the 6 j + ^ and 8j searches assuming SM 
branching ratios [42] for the W^* decay and is shown on Figure 9(a). The largest deviation 
1 — CLft = 0.0009 corresponding to 3.1a is reached at mH±=45 GeV and mA=44.9 GeV 
(in the middle of the very narrow lightest strip from mH±=40 GeV to 50 GeV at the 
{mii±,mA) diagonal). However, the mean background shift on the [mjj±,mA] plane amounts 
only to 1.1(7. 1 — CLh for the 4j -|- r channel is shown in Figure 9(b). The largest devi- 
ation 1 — CLj, = 0.013 corresponding to 2.2a" appears for low charged Higgs-boson masses 
(mH±=40 GeV, mA=21 GeV), reflecting the excess of events in the v^=189 GeV search. The 
mean background shift for this channel is O.Scr. Note that the results shown in Figures 9(a-b) 
are model-independent. 

When all channels are combined within 2HDM(I), a few hot spots with a significance 
above 2a survive. This is illustrated in Figures 9(c-d), where the combined results are plotted 
for tan /3=10 and 100. For tan (3=10, the largest excess 1 — CLb = 0.01 corresponding to 2.3a 
is found at mjj±=55 GeV and mA=35 GeV (just before switching from the bosonic to the 
fermionic channels). 

As mentioned previously, the H='=— s-AW='=* decay becomes dominant if the A boson is 
sufficiently light. The smaller tan /3 is, the smaller mA should be. This is clearly seen from the 
structure of the result on Figures 9(c-d): for tan/3=10, the bosonic decay becomes dominant 
at ruA ~ mH± — 18 GeV, while for tan/3=100, it dominates already at uia ~ mH±— 6 GeV. 

The model-independent limits on the charged Higgs-boson production cross section rela- 
tive to the 2HDM prediction are presented in Figure 10(a) for the H'*"H~^AW+*AW~* and in 
Figure 10(b) for H^H^^AW^*r^'T- searches, with the only assumption that W^* decays with 
SM branching ratios. The results combining all channels using 2HDM(I) branching ratios are 
shown on Figure 11 for different choices of tan/3. 
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OPAL: H"^ H" ^ AW* AW* OPAL: H"^ H"^AW*tv 




inH± (GeV) nijji (GeV) 



OPAL: H", 2HDM type I, tan|3=10 OPAL: H"^ H", 2HDM type I, tanp=100 




inH± (GeV) nijji (GeV) 



Figure 9: The conEdence 1 - CLh on the [mH±,mA] plane (a) for H+H~^AW+*AW"* 
combining the results of the 6j + £ and 8j searches, and (h) for H"'"H~^AW^*t^'t from the 
4j + T search. The combined results in 2HDM(I) for (c) tan P=10 and (d) tan P=100 are also 
shown. The significance values corresponding to the different shadings are shown by the bars 
at the right. 
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Figure 10: The 95% CL upper limits on the production cross section times relevant and A 
boson decay branching ratios relative to the 2HDM prediction on the [?n.jj± , itia] plane for the 
process (a) H+H-^AW+*AW-* and (b) H+H-^AW=^*rz^^. The plotted curves are isolines 
along which the limit is equal to the number indicated. The expected limits are given for (a) 
BR(H±^AW±*)^ • BR(^^bb)^ = 1 and (b) BR(H±^tz^^) • BR(H±^AW±*) • BR(yl^bb) = 
0.25 corresponding to the maximal value in 2HDM(I). 



For intermediate and large values of tan /3, an unexcluded region appears parallel to the 
[mii±,mA] diagonal both in the observed and expected limits in Figure 11 (see the full and 
dotted black curves corresponding to a cross-section ratio of 1.0). This can be understood 
by making two observations: Close to the diagonal, the eight- or six-jet structure of the 
H^^AW^* signal events becomes less pronounced and the final state turns four-jet-like. As 
the selection variables for the signal and the background become similar, the likelihood cut 
removes more signal events resulting in a drop in efficiency and therefore in sensitivity. This 
decreased sensitivity at the [mH±,"iA] diagonal is clearly visible in Figure 10(a). On the 
other hand, the rate of the fermionic events decreases in 2HDM(I) by the distance from the 
diagonal and by tan P, so the fermionic searches also lose their sensitivity. 

To further study the behavior of the unexcluded regions, 90%, 95% and 99% CL excluded 
areas are shown in Figure 12 for different choices of tan/3. At the 90% CL, the small unex- 
cluded islands at mA=12— 15 GeV disappear and the observed limit (as expected) is set at 
the transition region discussed in the previous paragraph. 

Due to the excess of events in the H"'"H~^AW"''*AW~* searches in the year 1999 data, 
the observed limit is lower than the expectation in all regions where the H^— >AW^* decay 
dominates. 

Finally our results are presented independent of tan/3 in Figure 13, and the limits on 
the charged Higgs-boson mass are summarized in Table 9. The absolute lower limit on the 
charged Higgs-boson mass is 56.5 GeV for tan/? < 100 and 12 GeV < mA < n-Ri) which 
should be compared to an expectation of 71.0 GeV. The observed limit is set by tan/?=3.5. 



24 



OPAL: H"^ H", 2HDM type I, tan (3=0.1 



_ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
L Ss^H'H-ya^HDM 




1 




1a): 


- — observed 












expected 


















0.6 


1.0 


2.0 J 






0.4 










0.2 


























40 



50 



60 



70 



80 



90 

■- (GeV) 



O 



90 
80 
70 
60 
50 
40 
30 
20 
10 



OPAL: H"^ H", 2HDM type I, tan (3=1 



40 



1 1 1 1 1 1 1 1 1 
L a,5(H*H-)/0 


2HDM 


' 1 


1 1 




- — observed 










expected 




0.6 


1.0 


12.0 ; 




0.2 loAJy 








: / 














50 



60 



70 



80 90 

mH±(GeV) 




mj,± (GeV) mj± (GeV) 



Figure 11: The 95% CL upper limits on the production cross section in 2HDM(I) relative to 
the theoretical prediction on the [mii±,mA] plane for different choices of tan/?; (a) 0.1, (b) 
1.0, (c) 10.0 and (d) 100.0. The plotted curves are isolines along which the limit is equal to 
the number indicated. 
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mj± (GeV) m^± (GeV) 




mj± (GeV) m^± (GeV) 



Figure 12: Excluded areas at 90%, 95% and 99% CL on the [m^±,mA] plane in 2HDM(I) for 
different clioices of tan 0: (a) 0.1, (b) 1.0, (c) 10.0 and (d) 100.0. 
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mH±(GeV) 

Figure 13: Excluded areas in 2HDM(I) on the [mH±,mA] plane independent of tan f3 at 95% 
CL. The weakest overall mass limit is defined by the tan (3=3.5 exclusion, which is also shown. 

where a small unexcluded island is present around mH±~60 GeV and mA~12 GeV (also 
shown on Figure 13). This unexcluded island is no longer present at 90% CL where the mass 
limit improves to 66.8 GeV. 

For toa>15 GeV, the tan /3-independent lower limit on the charged Higgs-boson mass at 
95% CL is 64.8 GeV with 71.0 GeV expected. Both the observed and expected limits are 
found in the transition region where the bosonic and fermionic channels have comparable 
sensitivities. The 6 GeV difference is due to the excess observed in the H"'"H~^AW''~*AW~* 
search. On the (mH±, toa) diagonal, the H^— >AW^* decay becomes kinematically suppressed, 
and the 2HDM(II) result for BR(H±^Ti/r)=0.65 of 81.6 GeV (see Figure 8 and Table 8) is 
reproduced. 

7 Summary 

A search is performed for the pair production of charged Higgs bosons in electron-positron 
collisions at LEP2, considering the decays H'^^tUt, qq and AW='=*. No signal is observed. 
The results are interpreted in the framework of Two-Higgs-Doublct Models. 

In 2HDM(II), required by the minimal super symmetric extension of the SM, charged Higgs 
bosons are excluded up to a mass of 76.6 GeV (with an expected hmit of 76.8 GeV) when 
BRiU^^TUr) + BR(H=^^qq) = 1 is assumed. See Figure 8 and Table 8 for BR(H=^^rz/^)- 
dependent limits. 

In the more exotic 2HDM(I), where fermionic decays can be suppressed and H^^AW^* 
can become dominant, a tan /3-independent lower mass limit of 56.5 GeV is observed for 
mA> 12 GeV (with an expected limit of 71.0 GeV) due to a non-signal-like excess observed 
at ^=192-202 GeV in the H+H"^AW+*AW~* search, discussed in Section 4. For mA> 
15 GeV, the observed limit improves to mH±> 64.8 GeV (with an expected limit of 71.0 
GeV). Figure 13 shows the excluded areas in the [mH±, wa] plane and Table 9 reports selected 
numerical results. 
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tan/3 


ruA 


limit on mH± (GeV) 






observed 


expected 


< 100 


12 GeV< m-A < ?n-H± 


56.5 (3.5) 


71.0 (1.0) 




mA=12 GeV 


56.5 (3.5) 


71.0 (1.0) 




mA=TOH± /2 


66.0 (3.5) 


73.5 (1.5) 




"lA > JTiHi — 10 GeV 


64.8 (100) 


71.7 (100) 




mA > JTijji — 5 GeV 


80.5 (100) 


78.0 (100) 


0.1 


12 GeV< m-A < "t-h* 


81.6 


81.1 




mA=12 GeV 


81.6 


81.1 




mp=mY{±/2 


81.8 


81.5 




ruA > ?TiH± — 10 GeV 


81.9 


81.6 




TTiA > Tnii± — 5 GeV 


81.9 


81.6 


1 


12 GeV< niA < ^h± 


66.4 


71.0 




mA=12 GeV 


66.4 


71.0 




mA=mii± /2 


78.3 


76.6 




mA > 'm}i± — 10 GeV 


81.9 


81.5 




mA > ?Wh± ~ 5 GeV 


81.9 


81.6 


10 


12 GeV< niA < 'tih* 


65.8 


73.9 




mA=12 GeV 


69.0 


82.5 




mA=?TiH± /2 


86.7 


89.8 




mA > JTT-Hi ~ 10 GeV 


81.5 


80.4 




niA > TiRi ~ 5 GeV 


81.9 


81.5 


100 


12 GeV< mA < 


64.8 


71.7 




mA=12 GeV 


69.3 


82.7 




mA=mH± /2 


87.1 


89.9 




niA > — 10 GeV 


64.8 


71.7 




TTiA > TiRi ~ 5 GeV 


80.5 


78.0 



Table 9: Lower mass limits for the charged Higgs boson in 2HDM(I). For the tan f3< 100 
results, the tan /? value at which the limit is set is indicated in parenthesis. For any tan /3 
value, an extrapolation of the exclusion hmits to mH± =mA gives the result quoted in Table 8 
for BR{Il^^TUr)=0.65. 
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Abstract 

A search is made for charged Higgs bosons predicted by Two-Higgs-Doublet extensions 
of the Standard Model (2HDM) using electron-positron collision data collected by the OPAL 
experiment at ^/s = 189 — 209 GeV, corresponding to an integrated luminosity of approxi- 
mately 600 pb~^. Charged Higgs bosons are assumed to be pair-produced and to decay into 
qq, TUr or AW^*. No signal is observed. Model-independent limits on the charged Higgs- 
boson production cross section are derived by combining these results with previous searches 
at lower energies. Excluded areas on the [mjji, BR(H^— ^rz^,-)] plane are presented assum- 
ing BR(H^^Ti/,-) + BR(H^^qq) = 1. Under the above assumption, motivated by general 
2HDM type II models, charged Higgs bosons are excluded up to a mass of 76.6 GeV at 95% 
confidence level, independent of the branching ratio BR(H^^r^'T-). A scan of the 2HDM type 
I model parameter space is performed and limits on the Higgs-boson masses mjji and ttt-a 
are presented for different choices of tan /?. 
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1 Introduction 

In the Standard Model (SM) [1], the electroweak symmetry is broken via the Higgs mecha- 
nism [2] generating the masses of elementary particles. This requires the introduction of a 
complex scalar Higgs-field doublet and implies the existence of a single neutral scalar parti- 
cle, the Higgs boson. While the SM accurately describes the interactions between elementary 
particles, it leaves several fundamental questions unanswered. Therefore, it is of great interest 
to study extended models. 

The minimal extension of the SM Higgs sector required, for example, by supersymmetric 
models contains two Higgs-field doublets [3] resulting in five Higgs bosons: two charged 
(H^) and three neutral. If CP-conservation is assumed, the three neutral Higgs bosons 
are CP-eigenstates: h and H are CP-even and A is CP-odd. Two-Higgs-Doublet Models 
(2HDMs) are classified according to the Higgs-fermion coupling structure. In type I models 
(2HDM(I)) [4], all quarks and leptons couple to the same Higgs doublet, while in type II 
models (2HDM(II)) [5], down-type fermions couple to the first Higgs doublet, and up-type 
fermions to the second. 

Charged Higgs bosons are expected to be pair-produced in the process e+e~— >^H+H~ at 
LEP, the reaction c~''c~— >^H^W^ having a much lower cross section [6]. In 2HDMs, the tree- 
level cross section [7] for pair production is completely determined by the charged Higgs-boson 
mass and known SM parameters. 

The H='= branching ratios are model-dependent. In most of the 2HDM(II) parameter space, 
charged Higgs bosons decay into the heaviest kinematically allowed fermions, namely rvr and 
quark pairs-*^. The situation changes in 2HDM(I), where the decay H='=— s-AW*"* can become 

^Throughout this paper charge conjugation is imphed. For simpUcity, the notation tUt stands for r'^i'r 
and T^i>T and qq for a quark and anti-quark of any flavor. 
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dominant if the ratio of the vacuum expectation values of the two Higgs-field doublets tan /3^1 
and the A boson is sufRcicntly light [8]. 

In this paper we search for charged Higgs bosons decaying into qq, tVt and AW^* using the 
data collected by the OPAL Collaboration in 1998—2000. The results are interpreted within 
general 2HDM(II) assuming BR(H='=— >tz/^) + BR(H='=— >-qq) = 1 for the branching ratios and 
in 2HDM(I) taking into account decays of charged Higgs bosons via AW^*, as well. Our 
result is not confined to qq={cs, cs} although that is the dominant hadronic decay channel 
in most of the parameter space. 

The previously published OPAL lower limit on the charged Higgs-boson mass, under the 
assumption of BR(H='=^rz/T-) + BR(H^^qq) = 1, is mjj± > 59.5 GeV at 95% confidence level 
(CL) using data collected at < 183 GeV [9, 10]. Lower bounds of 74.4 - 79.3 GeV have 
been reported by the other LEP collaborations [11-13] based on the full LEP2 data set. The 
DELPHI Collaboration also performed a search for H^— >AW^* decay and constrained the 
charged Higgs-boson mass in 2HDM(I) [12] to be mH±>76.7 GeV at 95%CL. 

2 Experimental considerations 

The OPAL detector is described in [14]. The events are reconstructed from charged-particle 
tracks and energy deposits {clusters) in the electromagnetic and hadron calorimeters. The 
tracks and clusters must pass a set of quality requirements similar to those used in previous 
OPAL Higgs-boson searches [15]. In calculating the total visible energies and momenta of 
events and individual jets, corrections are applied to prevent double-counting of energy in the 
case of tracks and associated clusters [15]. 

The data analyzed in this paper were collected in 1998—2000 at center-of-mass energies of 
189 — 209 GeV as given in Table 1. Due to different requirements on the operational state of 
the OPAL subdetectors, the integrated luminosity of about 600 pb~^ differs slightly among 
search channels. 



Year 


1998 


1999 


2000 


Ecm (GcV) 


186-190 


190-194 


194-198 


198-201 


201-203 


200-209 


<£cm> (GeV) 


188.6 


191.6 


195.5 


199.5 


201.9 


206.0 


(GeV) 


189 


192 


196 


200 


202 


206 


JCdt (pb"') (2t) 


183.5 


29.3 


76.4 


76.6 


45.5 


212.6 


jrdt (pb-^) (2j+r,4j) 


179.6 


29.3 


76.3 


75.9 


36.6 


217.4 


(pb-^) (8j, 6j+£,4j + r) 


175.0 


28.9 


74.8 


77.2 


36.1 


211.1 



Table 1: Data-taking year, center-of-mass energy bins, luminosity-weighted average centcr-of- 
mass energies, the energies of signal and background Monte Carlo simulations, and integrated 
luminosities of the data. The data correspond to total integrated luminosities of 623.9 pb~^ 
for the two-tau, 615.1 pb~^ for the two-jet plus tau and the four-jet channels and 603.1 ph"^ 
for the H='=— >AW='=* selections. 

In this paper the following final states are sought: 

• H+H~^r+^'T-''"~^'T {two-tau final state, 2r), 

• H+H~^qqrz/T- {two-jet plus tau final state, 2j + t), 

• H"*"H~^qqqq {four-jet final state, 4j), 

• H+H-^AW+*AW-*^bbqqbbqq {eight-jet final state, 8j), 
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• H+H-^AW+*AW"*^_bbqqbb£f£ {six-jet plus lepton final state, 6j +^), 

• H+H^ — >AW^*TZ/T-^bbqqT;/T- {four-jet plus tau final state, 4j + r). 

The signal detection efficiencies and accepted background cross sections are estimated 
using a variety of Monte Carlo samples. The HZHA generator [16] is used to simulate H+H" 
production at fixed values of the charged Higgs-boson mass in steps of 1 — 5 GeV from the 
kinematic limit down to 50 GcV for fermionic decays and 40 GeV for bosonic decays. 

The background processes arc simulated primarily by the following event generators: 
PYTHIA [17] and KK2F [18] (Z/7*^qq(7)), grc4f [19] (four-fermion processes, 4f), BH- 
WIDE [20] and TEEGG [21] (e+e-(7)), KORALZ [22] and KK2F (/x+/x-(7) and t+t~{j)), 
PHOJET [23], HERWIG [24], Vermaseren [25] (hadronic and leptonic two-photon processes). 

The generated partons, both for the signal and the SM Monte Carlo simulations, are 
hadronized using JETSET [17], with parameters described in [26]. For systematic stud- 
ies, cluster fragmentation implemented in HERWIG for the process Z/7*— >^qq(7) is used. 
The predictions of 4f processes are cross-checked using EXCALIBUR [27], KoralW [28] and 
KandY [29]. 

The obtained Monte Carlo samples are processed through a full simulation of the OPAL 
detector [30] . The event selection is described below. 

3 Search for four-fermion final states 

In most of the parameter space of 2HDM(II) and with a sufficiently heavy A boson in 
2HDM(I), the fermionic decays of the charged Higgs boson dominate and lead to four-fermion 
final states. The most important decay mode is typically H^^tUt, with the hadronic mode 
H^— >qq reaching about 40% branching ratio at maximum. 

The search for the fully leptonic final state H+H~— ^-t+z/t-t^P^ is described in [31]. The 
searches for the H+H~— ^^qqri^T- and the H+H~— >^qqqq events are optimized using Monte Carlo 
simulation of H^^cs decays. The sensitivities to other quark flavors arc similar and the 
possible differences are taken into account as systematic uncertainties. Therefore, our results 
are valid for any hadronic decay of the charged Higgs boson. 

Four-fermion final states originating from H+H" production would have very similar kine- 
matic properties to W^W~ production, which therefore constitutes an irreducible background 
to our searches, especially when mH± is close to m^i. To suppress this difficult SM back- 
ground, a mass-dependent likelihood selection (similar to the technique described in [32]) is 
introduced. For each charged Higgs-boson mass tested (mtest), a specific analysis optimized 
for a reference mass (jn^^.f) close to the hypothesized value is used. 

We have chosen a set of reference charged Higgs-boson masses at which signal samples are 
generated. Around these reference points, mass regions (labeled by mref ) are defined with the 
borders centered between the neighboring points. For each individual mass region, at each 
ccntcr-of-mass energy, we create a separate likelihood selection with the signal histograms 
built using events generated at micf. The background histograms are composed of the SM 
processes and are identical for all mass regions. 

When testing the hypothesis of a signal with mass mtestj the background and data rate 
and discriminant (i.e. the reconstructed Higgs-boson mass) distribution depend on the mass 
region to which m-tcst belongs. The signal quantities depend on the value of mtcst itself and are 
determined as follows. The signal rate and discriminant distribution are computed, with the 
likelihood selection optimized for mref, for three simulated signal samples with masses miow; 
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mrcf and mhigh- Here, miow and mhigh are the closest mass points to mref at which signal Monte 
Carlo samples are generated, with miow< nT-icfK "Tihigh. The signal rate and discriminant 
distribution for rritcst arc then calculated by linear interpolation from the quantities for miow 
and mref if "T'test< ^"rcf, or for rrircf and mhigh if "T'test> rnref- 

When building the likelihood function three event classes are considered: signal, four- 
fermion background and two-fermion background. The likelihood output gives the probability 
that a given event belongs to the signal rather than to one of the two background sources. 

3.1 The two-jet plus tau final state 

The analysis closely follows our published one at s/s = 183 GeV [10]. It proceeds in two 
steps. First, events consistent with the final state topology of an isolated tau lepton, a pair of 
hadronic jets and sizable missing energy are preselected and are then processed by a likelihood 
selection. The sensitivity of the likelihood selection is improved by building mass-dependent 
discriminant functions as explained above. 

Events are selected if their likelihood output {£) is greater than a cut value chosen to 
maximize the sensitivity of the selection at each simulated charged Higgs-boson mass. Apart 
from the neighborhood of the W+W~ peak, the optimal cut does not depend significantly on 
the simulated mass and is chosen to be jC> 0.85. Around the W"'"W~ peak, it is gradually 
reduced to 0.6 at the lowest. 

At the end of the selection, 331 events are selected in the data sample with 316.9 ± 

3.2 (stat.) ±38.4 (syst.) events expected from SM processes for a test mass of mH±=75 GeV. 
The sources of systematic uncertainties are discussed below. Four-fermion processes account 
for more than 99% of the SM background and result in a large peak in the reconstructed mass 
centered at the mass (with a second peak at the Z mass for test masses of mjji > 85 GeV). 
The signal detection efficiencies for the various LEP energies are between 25% and 53% for 
any charged Higgs-boson mass. 

The likelihood output and reconstructed di-jet mass distributions for simulated Higgs- 
boson masses of 60 GeV and 75 GeV are presented in Figures l(a-d). The reconstructed 
Higgs-boson mass resolution is 2.0 — 2.5 GeV [10]. Figure 2(a) gives the mass dependence of 
the expected number of background and signal events and compares them to the observed 
number of events at each test mass. 

The systematic uncertainties are estimated for several choices of the charged Higgs-boson 
mass from 50 GeV to 90 GeV at ccnter-of-mass energies of ■y/i=189 GeV, 200 GeV and 
206 GeV to cover the full LEP2 range. The following sources of uncertainties are considered: 
limited number of generated Monte Carlo events, statistical and systematic uncertainty on 
the luminosity measurement, modeling of kinematic variables in the pre- and likelihood se- 
lections, tau lepton identification, dependence of the signal detection efficiency on final-state 
quark flavor, signal selection efficiency interpolation between generated Monte Carlo points, 
background hadronization model, and four-fermion background model. The contributions 
from the different sources are summarized in Table 2. 

In the limit calculation, the efficiency and background estimates of the 2] + t channel 
are reduced by 0.8—1.7% (depending on the center-of-mass energy) in order to account for 
accidental vetoes due to accelerator-related backgrounds in the forward detectors. 
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OPAL: H*H"^qqTV, 189-209 GeV 
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Figure 1: Likelihood output and reconstructed di-jet mass distributions for the (a-d) 2 j + r 
and (e-h) 4j channels. The distributions are summed up for all center-of-mass energies and 
correspond to 60 GeV and 75 GeV simulated charged Higgs-boson masses. All Monte Carlo 
distributions are normalized to the integrated luminosity of the data. When plotting the 
likelihood output, the signal expectation is scaled up by a factor of 10 for better visibility. A 
hadronic branching ratio of 0.5 is assumed for the 2 j + r signal, and 1.0 for the 4j signal. The 
reconstructed mass distributions are shown after the likelihood selection. 
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OPAL: H*H qqxv, 189-209 GeV OPAL: H*H qqqq, 189-209 GeV 




50 60 70 80 90 100 50 60 70 80 90 



m„± (GeV) m„± (GeV) 

Figure 2: The number of observed data, expected background and signal events for the (a) 
2 j + r and (b) 4j channels. The numbers are summed up for all center-of-mass energies and 
shown as a function of the reference charged Higgs-boson mass. A hadronic branching ratio 
of 0.5 is assumed for the 2j + r signal, and 1.0 for the 4j signal. Each bin corresponds to a 
different likelihood selection optimized for the mass at which the dot is centered. Since the 
same background simulations are used to form the reference histograms and the same data 
enter the selection, the neighboring points are strongly correlated. 



Source 


2j+r 


4j 


signal 


background 


signal 


background 


MC statistics 


3.1-4.6 


1.4-4.3 


1.6-2.4 


0.9-1.9 


luminosity 


0.3 


0.3 


0.3 


0.3 


preselection 


1.5-4.7 


1.8-7.6 


0.3-1.1 


0.5-2.2 


likelihood selection 


0.9-6.5 


5.8-22.7 


0.7-2.4 


2.1-7.5 


tau identification 


3.0 


3.0 


N.A. 


N.A. 


quark flavor 


2.7-3.8 


N.A. 


1.2-6.4 


N.A. 


interpolation 


0.2-0.4 


N.A. 


0.7-3.7 


N.A. 


hadronization model 


N.A. 


1.0-2.7 


N.A. 


0.7-4.1 


4f background model 


N.A. 


0.3-3.3 


N.A. 


1.7-3.7 



Table 2: Relative systematic uncertainties on the expected background and signal rates for 
the 2j + r and 4j Gnal states. The numbers are given in % and depend on the centre-of-mass 
energy and the reference charged Higgs-boson mass. N.A. stands for not applicable. 
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3.2 The four-jet final state 

The event selection follows our published analysis at i/i=183 GeV [10]: first, well-separated 
four-jet events with large visible energy are preselected; then a set of variables is combined 
using a likelihood technique. To improve the discriminating power of the likelihood selec- 
tion, a new reference variable is introduced: the logarithm of the matrix element probability 
for W+W~ production averaged over all possible jet-parton assignments computed by EX- 
CALIBUR [27]. Moreover, we introduce mass-dependent likelihood functions as explained 
above. As the optimal cut value on the likelihood output is not that sensitive to the charged 
Higgs-boson mass in this search channel, we use the condition L > 0.45 at all center-of-mass 
energies and for all test masses. 

There is a good agreement between the observed data and the SM Monte Carlo expec- 
tations at all stages of the selection. After all cuts, 1100 events are selected in the data, 
while 1117.8 ± 5.9 (stat.) ±74.4 (syst.) events are expected from SM processes for a test 
mass of mH±=75 GeV. The four-fermion processes account for about 90% of the expected 
background and result in a large peak centered at the mass and a smaller one at the Z 
boson mass. The signal detection efficiencies are between 41% and 59% for any test mass and 
center-of-mass energy. 

Typical likelihood output and reconstructed di-jet mass distributions of the selected events 
together with the SM background expectation and signal shapes for simulated charged Higgs- 
boson masses of 60 GeV and 75 GeV are plotted in Figures l(e-h). The Higgs-boson mass 
can be reconstructed with a resolution of 1 — 1.5 GeV [10]. Figure 2(b) shows the mass 
dependence of the expected number of background and signal events and compares them to 
the observed number of events at each test mass. Systematic uncertainties are estimated in 
the same manner as for the 2j +t search and are given in Table 2. 

4 Search for AW+*AW~* events 

In a large part of the 2HDM(I) parameter space, the branching ratio of H^ — ^ AW^* dom- 
inates. The possible decay modes of the A boson and the W^* lead to many possible 
H+H~^AW"'"*AW~* event topologies. Above mA~12 GeV, the A boson decays predomi- 
nantly into a bb pair, and thus its detection is based on b-flavor identification. Two pos- 
sibilities, covering 90% of the decays of two W='=*, arc considered: quark pairs from both 
W^* bosons or a quark pair from one and a Icptonic final state from the other. The event 
topologies are therefore "eight jets" or "six jets and a lepton with missing energy", with four 
jets containing b-fiavor in both cases. 

The background comes from several Standard Model processes. ZZ and W+W" produc- 
tion can result in multi-jet events. While ZZ events can contain true b-flavored jets, W^W~ 
events are selected as candidates when c-flavored jets fake b-jets. Radiative QCD corrections 
to e"'"e~-^qq also give a significant contribution to the expected background. 

Due to the complexity of the eight-parton final state, it is more efficient to use general event 
properties and variables designed specifically to discriminate against the main background 
than a full reconstruction of the event. As a consequence, no attempt is made to reconstruct 
the charged Higgs-boson mass. 

The analysis proceeds in two steps. First a preselection is applied to select b-tagged multi- 
jet events compatible with the signal hypothesis. Then a likelihood selection (with three event 
classes: signal, four-fermion background and two-fermion background) is applied. 
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The preselection of multi-jet events uses the same variables as the search for the hadronic 
final state in [10] with optimized cut positions. However, it introduces a very powerful new 
criterion, especially against the W+W" background, on a combined b-tagging variable (^Sevt) 
requiring the consistency of the event with the presence of b-quark jets. 

The neural network method used for b-tagging in the OPAL SM Higgs-boson search [15] 
is used to calculate on a jet-by-jet basis the discriminating variables f^^^ and /^^g/^- These 
are constructed for each jet i as the ratios of probabilities for the jet to be c- or uds-like 
versus the probability to be b-like. The inputs to the neural network include information 
about the presence of secondary vertices in a jet, the jet shape, and the presence of leptons 
with large transverse momentum. The Monte Carlo description of the neural network output 
was checked with LEPl data with a jet energy of about 46 GeV. The main background in 
this search at LEP2 comes from four-fermion processes, in which the mean jet energy is about 
50 GeV, very close to the LEPl jet energy; therefore, an adequate modeling of the data is 
expected with an event reconstruction assuming four jets. 

The AW+*AW~* signal topology depends on the Higgs-boson masses. At mA~12 GeV 
or mA«mH±, the available energy in the A or W^* system is too low to form two clean, 
coUimated jets. At high mH± , the boost of the A and W='=* bosons is small in the laboratory 
frame and the original eight partons cannot be identified. At low mjji, the A and W^* 
bosons might have a boost, but it is still not possible to resolve correctly the two partons 
from their decay. Prom these considerations, one can conclude that it is not useful to require 
eight (or even six) jets in the event, as these jets will not correspond to the original partons. 
Consequently, to get the best possible modeling of the data, four jets are reconstructed with 
the Durham jet-finding algorithm [33] before the b-tagger is run. 

The flavor-discriminating variables are combined for the four reconstructed jets by 

The index i runs over the reconstructed jets (f = 1, ...4) and the parameters a and (3 are 
numerical coefficients whose optimal values depend on the flavor composition of the signal 
and background final states. However, since the expected sensitivity of the search is only 
slightly dependent on the values of a and (3, they are fixed at a = 0.1 and /3 = 0.7. Events 
are retained if Bevt > 0.4. 

The preselections of the two event topologies (8j and 6j + i) arc very similar. However, in 
the 6i+£ channel, no kinematic fit is made to the W+W^-^qqqq hypothesis and, therefore, 
no cuts are made on the fit probabilities. No lepton identification is applied; instead the 
search is based on indirect detection of the associated neutrino by measuring the missing 
energy. 

After the preselection the observed data show an excess over the predicted Monte Carlo 
background. This can partly be explained by the apparent difference between the gluon 
splitting rate into cc and bb pairs in the data and in the background Monte Carlo simulation. 
The measured rates are g^c = 3.2 ± 0.21 ± 0.38% [34] and ^^B = 0-307 ± 0.053 ± 0.097% [35] 
from the LEPl OPAL data. The gluon splitting rates in the Monte Carlo are extracted from 
e+e-^ZZ^£+^-qq events and are found to be g^'^ = 1.33 ± 0.06% and g^ = 0.116 ± 
0.0167%, averaged over all center-of-mass energies. This mismodeling can be compensated 
by reweighting the SM Monte Carlo events with gluon splitting to heavy quarks by universal 
reweighting factors [36] and at the same time deweighting the non-split events to keep the 
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total numbers of W+W", ZZ and two-fermion background events fixed at generator level. 
The reweighting factor is 2.41 for g^cc and 2.65 for g^bb. This correction results in a 
background enhancement factor of 1.08 to 1.1 after the preselection, depending on the search 
channel and the center-of-mass energy, but it does not affect the shape of the background 
distributions. The numbers of preselected events after the reweighting are given in Table 3. 



LEP energy 


preselection 


preselection 


exclusive 


exclusive 


overlap 


(year) 


8j 


Gj + (. 


8j 


Gj + /' 




189 GeV data 


238 


358 


3 


24 


5 


(1998) background 


231.2±2.9 


342.2±3.6 


2.1±0.3 


24.4±1.0 


6.3±0.5 


192 - 202 GeV data 


297 


310 


16 


16 


17 


(1999) background 


270.4±2.9 


285.0±3.0 


13.3±0.7 


10.4±0.6 


13.4±0.7 


200 - 209 GeV data 


265 


281 


9 


8 


15 


(2000) background 


252.5±3.7 


270.5±5.0 


13.0±0.9 


9.3±0.8 


12.9±0.9 



Table 3: Observed data and expected SM background events for each year in the AW"'"*AW~* 
searches. The 8j and 6 j + £ event samples after the preselection step (2nd and 3rd columns) 
are highly overlapping. After the likelihood selection, the overlapping events are removed 
from the 8j and 6j + ^ samples and form a separate search channel (last three columns). The 
uncertainty on the background prediction due to the limited number of simulated events is 
given. The Monte Carlo reweighting to the measured gluon splitting rates is included. 

As a final selection, likelihood functions are built to identify signal events. The reference 
distributions depend on the LEP energy, but they are constructed to be independent of the 
considered (mH± ) "^a) combination. To this end, we form the signal reference distributions 
by averaging all simulated H+H" samples in the {m}i±,mA) mass range of interest. 

Since the selections at ^/s = 192 — 209 GeV are aimed at charged Higgs-boson masses 
around the expected sensitivity reach of about 80—90 GeV, all masses up to the kinematic 
limit are included. On the other hand, at v^=189 GeV only charged Higgs-boson masses 
up to 50 GeV are included since the selections at this energy are optimized to reach down 
to as low as a charged Higgs-boson mass of 40 GeV where the LEPl exclusion limit lies. 
The input variables for the 8j final state are: the Durham jet-resolution parameters^ logio 2/34 
and logioy56) the oblateness [37] event shape variable, the opening angle of the widest jet 
defined by the size of the cone containing 68% of the total jet energy, the charge-signed 
cosine of the production angle in the W^W"— ^qqqq hypothesis, and the b-tagging variable 
Bcvt- At -^5=189 GeV, log^gyss, logioy45, logioye?) and the maximum jet energy are also 
used. Moreover, the sphericity [38] event shape variable has more discriminating power and 
thus replaces oblateness. Although the yij variables are somewhat correlated, they contain 
additional information: their differences reflect the kinematics of the initial partons. 

The input variables for the 6] + i selection are: log 2/34, log^Qj/sg, the oblateness, the 
missing energy of the event, and i3evt- At y^=189 GeV, logj^gysa, the maximum jet energy 
and the sphericity are also included. 

Events are selected if they pass a lower cut on the likelihood output. This final cut 
docs not rcmow (especially in the year 1999 data) all the excess events observed after the 

"^Throughout this paper yij denotes the parameter of the Durham jet finder at which the event classification 
changes from i-jet to j-jet, where j = i + 1. 
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preselection. The distributions of the critical selection variables, 2/34, y^Q and i3evt, are plotted 
on Figure 3 both in background-enriched data samples and after the preselection. To prepare 
these background-enriched data samples, the preselection cuts on 1/34 and B^yt dropped, 
except for the study of the yse variable where we keep the cut on 1/34 in order to select multi-jet 
events. The resulting samples are completely dominated by background processes. We find 
systematic differences between the data samples and the equivalent background Monte Carlo 
simulations at both stages. The observed excess cannot be attributed to a Higgs-boson signal 
whose contribution is expected to be at most 0.5% in the background-enriched data samples. 
The interpretation of the excess in terms of a systematic uncertainty on the background 
prediction will be discussed later. 

The non-signal nature of the excess is also demonstrated by Figure 4, which shows the 
distributions of the likelihood output. The positions of the likelihood cut are indicated by 
vertical lines. 

Because the 8j and the 6] + i selections have several discriminating variables in common, 
a major overlap between these selections exists. To assure that every event is counted only 
once, the two samples are redistributed into three: (i) events exclusively classified as 8j 
candidates, (m) events exclusively classified as 6j -I- ^ candidates and (in) events accepted 
by both selections. If an event falls into class (Hi), the larger likelihood output of the two 
selections is kept for further processing. The final results using the above classification are 
quoted in Table 3. 

This modified channel definition not only removes the overlap but also increases the ef- 
ficiency for detecting signal events by considering the cross-channel efficiencies (e.g. the 
efficiency to select H"'"H^^bbqqbbqq signal by the exclusive 6] + i selection can be as high 
as 18%, though it is typically only a few %). The efficiencies are determined independently 
for all simulated (mH±,wiA) combinations and interpolated to arbitrary (mH±,mA) by two- 
dimensional spline interpolation. The behavior of the selection efficiencies depends strongly 
on the targeted charged Higgs-boson mass range and also varies with the mass difference 
Am = mH± — ruA- In most cases the overlap channel has the highest efficiency. At y^=189 
GeV and mH±=45 GeV, it reaches 32% for the bbqqbb^f^ and 44% for the H+H~— ^^bbqqbbqq 
signal close to the mjj±=mA diagonal. At ■y/s=206 GeV and mH±=90 GeV, the overlap effi- 
ciency can be as high as 62% for the bbqqbb^f^ and 71% for the H+H~-^bbqqbbqq signal. 
The exclusive 6j -I- ^ selection has efficiencies typically below 20—30%, while the exclusive 8j 
selection below 10—15%. Table 4 gives the selection efficiencies at selected {m^±,mA) points. 

The composition of the background depends on the targeted Higgs-boson mass region. 
In the low-mass selection (■y/i=189 GeV) that is optimized for mjj±=40— 50 GeV, the Higgs 
bosons are boosted and therefore the final state is two-jet-like with the largest background 
contribution coming from two-fermion processes: they account for 52% in the exclusive 8j, 
80% in the exclusive + i and 76% in the overlap channel. On the other hand, in the high- 
mass analysis (-^5 = 192 — 209 GeV) the four-fermion fraction is dominant: it is 69% in the 
8j, 56% in the 6 j + ^ and 70% in the overlap channel. 

Although our studies show that the observed excess does not lie preferentially in the 
phase space region where the Higgs-boson signal is expected and that it probably originates 
from deficiencies of the Monte Carlo description of the jet-resolution parameters yij and the 
b-tagging variable B^vt, no further correction is applied to the estimated background in the 
background subtraction procedure of the statistical analysis to calculate exclusion limits. This 
leads to conservative limits on the production cross section. Therefore, no systematic uncer- 
tainty is assigned to the modeling of yu and Bevt at the preselection level. The systematic 
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OPAL: H^H^bbqqbbqq, 192-209 GeV 
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Figure 3: Most important selection variables: (a-b) logxoy34, (c-d) log^o 2/56 and (e-f) Bevt in 
the 8j channel at ^/s = 192 — 209 GeV. The distributions are shown (left) in a background- 
enriched data sample (see text for explanation) and (right) after the full preselection. To 
form the signal histograms, the Monte Carlo distributions are averaged for all simulated 
(mjj±,mA) mass combinations in the mass range of interest. The Monte Carlo reweighting 
to the measured gluon splitting rates is included. The expectations from SM processes are 
normalized to the data luminosity. The preselection cuts on 2/34 and Be^t are indicated by 
vertical lines. 13 




OPAL: ff'H^ bbqqbblv, 189 GeV OPAL: H^ff^ bbqqbblv, 192-209 GeV 




Likelihood Likelitiood 



Figure 4: Likelihood output distributions for the (a-b) 8j and (c-d) 6] + i channels at 
y/s=189 GeV and 192 — 209 GeV. To form the signal histograms, the Monte Carlo dis- 
tributions are averaged for all simulated (mH±,mA) mass combinations in the mass range of 
interest. The Monte Carlo reweighting to the measured gluon splitting rates is included. The 
expectations from SM processes are normalized to the data luminosity. The lower cuts on the 
likelihood output are indicated by vertical lines. 
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(mH±,mA) 


(GcV, GcV) 






(45,30) 


(80,50) 


(45,30) 


(90,60) 






^=189 GeV 


^=206 GcV 


bbqqbbqq 


8j 


4.6 


1.0 


9.3 


12.4 




overlap 


41.0 


2.8 


14.9 


69.6 




6j+£ 


17.0 


3.6 


4.1 


3.1 




total 


62.6 


7.4 


28.3 


85.0 


bbqqbb^t'^ 


6j+£ 


28.2 


6.0 


11.6 


7.0 




overlap 


31.8 


3.6 


14.2 


62.2 




8j 


1.8 


0.1 


2.2 


6.1 




total 


61.8 


9.7 


28.0 


75.3 


bbqqri^T- 


4j + r 


68.0 


0.0 


12.3 


11.1 



Table 4: Signal selection efRciencies in percent for the H^^AW^* final states in the different 
search channels at ^/s=189 and 206 GeV at representative (mH±,?TiA) points. 

uncertainties related to the other preselection variables, estimated from background-enriched 
data samples, are taken into account. The effect of mismodeling the shapes of the reference 
distributions in the likelihood selection is estimated for all variables, including and Bevt 
and is accounted for in the statistical analysis. 

Systematic uncertainties arise also due to the gluon splitting correction: the experimental 
uncertainty on the gluon splitting rate translates into uncertainties on the total background 
rates. Moreover, the weighted background counting introduces an uncertainty due to the 
Monte Carlo statistics of the g^cc and bb events. 

In summary, the following sources of systematic uncertainties are considered: limited 
number of simulated signal and background events, modeling of the preselection variables 
(other than y^j and Bg^t), modeling of the shapes of the reference distributions in the likelihood 
selection and the gluon splitting correction. Uncertainties below the 1% level are neglected. 
The different contributions are summarized in Table 5. 

5 Search for AW^*TUr events 

In some parts of the 2HDM(I) parameter space, both the fermionic H^^rz/,- and the bosonic 
11=*=— >AW^* decay modes contribute. To cover this transition region parallel to the (mH± , mA) 
diagonal, a search for the final state H+H~— >^AW='=*tz/t- is performed. The transition region 
is wide for small tan/3 and narrow for large tan/3; therefore, this analysis is more relevant for 
lower values of tan/3. 

Only the hadronic decays of W^* and the decay A^bb are considered. Thus the events 
contain a tau lepton, four jets (two of which are b-flavored) and missing energy. Separating 
the signal from the W^W~ background becomes difficult close to m}i±=m-\^±. 

The preselection is designed to identify hadronic events containing a tau lepton plus 
significant missing energy and transverse momentum from the undetected neutrino. In most 
cases it is not practical to reconstruct the four jets originating from the AW='=* system. Instead, 
to suppress the main background from semi-leptonic W"'"W~ events, we remove the decay 
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Source 


exclusive 8j 


exclusive 6j + ^ 


overlap 


signal 


background 


signal 


background 


signal 


background 


MC statistics 


>15/>8.5 


13.2/6.7 


>5.7/>8.4 


4.0/8.3 


>4.0/>2.8 


7.9/7.1 


preselection 


1.0 


1.0 


1.0 


2.0 


1.0 


2.0/1.5 


L selection 
















4.0/1.8 


6.0/6.2 


2.2/2.6 


6.0/8.0 


2.8/1.5 


5.5/4.9 


b-tag 


0.0/1.8 


4.7/7.0 


2.9/1.4 


4.4/7.1 


1.0/1.3 


4.0/5.1 


other 


1.8/0.7 


3.9/3.2 


1.0/1.6 


3.7/3.5 


1.2/0.7 


3.4/2.4 


gluon splitting 














g— >cc, exp. 


N.A. 


0.6/2.5 


N.A. 


1.6 


N.A. 


1.4/2.3 


g^cc, MC 


N.A. 


0.2/0.8 


N.A. 


0.6 


N.A. 


0.5/0.8 


g^bb, exp. 


N.A. 


1.4/4.2 


N.A. 


3.8/4.3 


N.A. 


5.5/5.4 


g^bb, MC 


N.A. 


0.5/1.7 


N.A. 


1.5/1.7 


N.A. 


2.2 


gluon splitting 
correction factor 


N.A. 


1.05/1.18 


N.A. 


1.13/1.15 


N.A. 


1.15/1.22 



Table 5: Relative systematic uncertainties in percent for the AW'''*AW~* searches. Where 
two values are given separated by a "/", the first belongs to the 189 GeV selection and the 
second to the 192 — 209 GeV selections. For the signal, the uncertainties due to the limited 
Monte Carlo statistics are calculated by binomial statistics for a sample size of 500 events 
and they also depend, via the selection efficiency, on the assumed Higgs-boson masses. N.A. 
stands for not applicable. The multiplicative gluon splitting correction factors, used to obtain 
the background-rate estimates as explained in the text, are given in the last line. 

products of the tau candidate and force the remaining hadronic system into two jets by the 
Durham algorithm. The requirements are then based on the preselection of Section 3.1 with 
additional preselection cuts on the effective center-of-mass energy, logj^Q j/12 and logj^g 2/23 of 
the hadronic system, and the charge-signed production angle. 

The Hkehhood selection uses seven variables: the momentum of the tau candidate, the 
cosine of the angle between the tau momentum and the nearest jet, log^^Q yi2 of the hadronic 
system, the cosine of the angle between the two hadronic jets, the charge-signed cosine of 
the W='= production angle, the invariant mass of the hadronic system, and the b-tagging 
variable Bc^t- Here, Bcvt is defined using the two jets of the hadronic system using Eq. (1) 
of Section 4, with i = 1,2 and a = f3 = I. To form the signal reference distributions, all 
simulated H'^H" samples in the (mH± , wa) mass range of interest are summed up. Since the 
search at i/s = 192—209 GeV targets intermediate charged Higgs-boson masses (60—80 GeV), 
all masses up to the kinematic limit are included. At v^=189 GeV, only charged Higgs-boson 
masses up to 50 GeV are included since the selection is optimized for low charged Higgs-boson 
masses (40-50 GeV). 

The likelihood output distributions are shown in Figure 5. There is an overall agreement 
between data and background distributions, apart from a small discrepancy at 189 GeV. 
Events are selected if their likelihood output is larger than 0.9. In total, 15 data events survive 
the selection at ^/s = 192 - 209 GeV, to be compared with 14.8 ± 0.6 (stat.) ±1.9 (syst.) 
events expected from background sources. At y^=189 GeV, where the selection is optimized 
for low Higgs-boson masses, 13 data events are selected with 6.1 it 0.5 (stat.) ±1.3 (syst.) 
events expected. The contribution of four-fermion events, predominantly from semi-leptonic 
W+W" production, amounts to 33% at ^^=189 GeV and to 90% at ^/s = 192 - 209 GeV. 
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Likelihood Likeliliood 

Figure 5: Likelihood output distribution for the 4j + r channel at (a) y/s=189 GeV and 
(b) ^/s=192 — 209 GeV. To form the signal histograms, the Monte Carlo distributions are 
averaged for all simulated {m^± , mA) mass combinations in the mass range of interest. The 
expectations from SM processes are normalized to the data luminosity. The lower cuts on the 
likelihood output are indicated by vertical line. 

At ^/s = 192 — 209 GeV, the signal selection efficiency starts at about 5% at mH±=40 GeV, 
reaches its maximum of about 40% (depending on the mass difference Am = — m.^) 
at mH±=60 GeV, then decreases to 12% at mH±=90 GeV. In the low- mass selection at 
■^5=189 GeV, the efficiency depends strongly on the mass difference: at mH±=40 GeV, it is 
22% for Am = 2 GeV and 60% for Am = 10 GeV. The selection efficiency approaches its 
maximum at mH±=50 GeV (70% for Am = 10 GeV) and then drops to zero at mjj±=80 GeV. 
Table 4 gives selection efficiencies at representative {mii±,mA) points. 

The systematic uncertainties due to the modeling of selection variables are evaluated with 
the method developed for the AW+*AW~* channels and summarized in Table 6. 



Source 


4j +r 




signal 


background 


MC statistics 


>2.7/>4.5 


8.2/7.0 


preselection: tau ID 


0.0 


2.3/5.0 


other 


0.0/1.0 


9.5/7.9 


likelihood selection: b-tag 


0.3/1.4 


2.4/3.0 


other 


3.2/2.1 


18.4/11.2 



Table 6: Systematic uncertainties in percent for the 4j + r channel. Where two values are 
given separated by a "/", the first one belongs to the 189 GeV selection and the second to the 
192 — 209 GeV selections. For the signal, the uncertainties due to the limited Monte Carlo 
statistics are calculated by binomial statistics for a sample size of 500 events and they also 
depend, via the selection efRciency, on the assumed Higgs-boson masses. 
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6 Interpretation 



None of the searches has revealed a signal-hke excess over the SM expectation. The re- 
sults presented here and those published previously [10, 31] by the OPAL Collaboration are 
combined using the method of [39] to study the compatibility of the observed events with 
"background-only" and "signal plus background" hypotheses and to derive limits on charged 
Higgs-boson production. The statistical analysis is based on weighted event counting, with 
the weights computed from physical observables, also called discriminating variables of the 
candidate events (see Table 7). Systematic uncertainties with correlations are taken into ac- 
count in the confidence level (p- value) calculations. To improve the sensitivity of the analysis, 
they are also incorporated into the weight definition [39]. 



Channel 


(GeV) 


Discriminant 


2r 


183 


simple event counting 


2r 


189-209 


likelihood output 


2j + r 


183-209 


reconstructed di-jet mass 


4j 


183-209 


reconstructed di-jet mass 


8j 


189 


simple event counting 


8j 


192-209 


likelihood output 


6j + £ 


189 


simple event counting 


6j + £ 


192-209 


likelihood output 


4j+T 


189-209 


simple event counting 



Tabic 7: Discriminating variables entering the statistical analysis for each search topology. 
Previously published results are also included. 

The results are interpreted in two different scenarios: in the traditional, super symmetry- 
favored 2HDM(II) (assuming that there are no new additional light particles other than the 
Higgs bosons) and in the 2HDM(I) where under certain conditions fermionic couplings are 
suppressed. 

First, we calculate 1 — CLf,, the confidence [39] under the background-only hypothesis, and 
then proceed to calculate limits on the charged Higgs-boson production cross section in the 
signal -I- background hypothesis. These results are used to provide exclusions in the model 
parameter space, and in particular, on the charged Higgs-boson mass. 

2HDM Type II 

First a general 2HDM(II) is considered, where BR(H^^rz/^) -|- BR(H^^qq) = 1. This model 
was thoroughly studied at LEP. It is realized in super symmetric extensions of the SM if no 
new additional light particles other than the Higgs bosons are present. As our previously 
published mass limit in such a model is mH± > 59.5 GeV [10], only charged Higgs-boson 
masses above 50 GcV are tested. Cross-section limits for lower masses can be found in [9]. In 
this model, the results of the 2r, 2] + t and 4j searches enter the statistical combination. 

The confidence 1 — CLh is plotted for each channel separately in Figure 6(a) and com- 
bined in Figure 6(b). Note that 1 — CLj, < 0.5 translates to negative values of sigma (as 
indicated by the dual y-axis scales in Figure 6(a)) and indicates an excess of events. The 
largest deviation, a 2.8a excess observed in the 2j + t channel at mH±=88 GeV, comes 
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from a known deficiency [40] of the Monte Carlo simulation of isolated tracks from the frag- 
mentation and hadronization process. When the three final states are combined assuming 
BR(H^— >tz/t-) + BR(H^— >qq) = 1, no deviation reaches the 2a level. 




Figure 6: The observed confidence levels for the background interpretation of the data, 1 — 
CLb, (a) for the three different final states as a function of the charged Higgs-boson mass, and 
(b) for the combined result in 2HDM(II) assuming BR(H=*=^TZ^r) + BR(H±^qq) = 1 on the 
[m,jj±, BR(H^^rz^T-)] plane. The significance values corresponding to the different shadings 
are shown by the bar at the right. 

The results are used to set upper bounds on the charged Higgs-boson pair production 
cross section relative to the 2HDM prediction as calculated by HZHA. The limits obtained 
are shown for each channel separately in Figures 7(a-c) and combined in Figure 7(d). The 
combined results are shown by "isolines" along which iT95(H+H~)/(T2hdm, the ratio of the 
limit on the production cross section and the 2HDM cross-section prediction, is equal to the 
number indicated next to the curves. 

Excluded areas on the [mjj± , BR(H^^r^'T-)] plane are presented for each channel sepa- 
rately in Figure 8(a) and combined in Figure 8(b). The expected mass limit from simulated 
background experiments, assuming no signal, is also shown. For the combined results, the 
90% and 99% CL contours are also given. Charged Higgs bosons are excluded up to a mass 
of 76.6 GeV at 95% CL, independent of BR(H±^rz/^). Lower mass limits for different values 
of BR(H^— >Ti/^) are presented in Table 8. 
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OPAL: H"^ H" ^ TVTV OPAL: H"^ H" qqxv 




mjj± (GeV) mjj± (GeV) 

Figure 7: Observed and expected 95% CL upper limits on the H+H^ production cross section 
times the relevant decay branching ratios relative to the theoretical prediction for the (a) 
rvT-rVr (b) qqTiy-r and (c) qqqq channels. The horizontal lines indicate the maximum possible 
branching ratios for a given channel. In (b), BKqqri/ = 2 • BR(H^^t^'t-) • BR(H^— >qq). 
(d) Upper limits on the production cross section relative to the 2HDM prediction on the 
[mH±, BR(H=''^rz^T-)] plane in 2HDM(n) assuming BR(H^^TUr) + BR(H±^qq) = 1. The 
plotted curves are isolines along which the limit is equal to the number indicated. 
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Figure 8: Observed and expected excluded areas at 95% CL on the [mH± , BR(H^— ^ri/,-)] 
plane (a) for each search channel separately and (h) combined in 2HDM(n) assuming 
BR(H=^^rz^^) + BR(H=^^qq) = 1. For the combined result, the 90% and 99% CL observed 
limits are also shown. See Table 8 for numerical values of the combined limit. 





Lower mass 


limit (GeV) 


Observed 


Expected 





76.9 


77.9 


0.5 


79.2 


78.0 


0.65 


82.0 


81.7 


1 


91.2 


89.2 


any 


76.6 (0.15) 


76.8 (0.27) 



Table 8: Observed and expected lower limits at 95% CL on the mass of the charged Higgs 
boson in 2HDM(II) assuming liI{(H^^Tiyr) + BR(H^^qq) = 1. For the results independent 
of the branching ratio (last line), the BR(H^— ^rz^,-) vaiue at which the limit is set, is given in 
parenthesis. 
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2HDM Type I 

We present here for the first time an interpretation of the OPAL charged Higgs-boson searches 
in an alternative theoretical scenario, a 2HDM(I). The novel feature of this model with respect 
to the more frequently studied 2HDM(II) is that the fermionic decays of the charged Higgs 
boson can be suppressed. If the A boson is light, the H^^AW^* decay may play a crucial 
role. 

The charged Higgs-boson sector in these models is described by three parameters: mjji, 
ruA and tan/3. To test this scenario, the charged Higgs-boson decay branching ratios are 
calculated by the program of Akeroyd et al. [8], and the model parameters are scanned in 
the range: 40 GeV < mH± < 94 GeV, 12 GeV < niA < mH±, 0.1 < tan/3 < 100. Charged 
Higgs-boson pair production is excluded below 40 GeV by the measurement of the Z boson 
width [41]. As the A boson detection is based on the identification of b-quark jets, no limits 
are derived for mA< 2m\). Below tan/3=0.1, BR(H^^AW^*) vanishes and the limit is no 
longer sensitive to ruA- 

Both the fermionic (2r, 2j + t and 4j) and the bosonic (4j +t, 6i + £ and 8j) final states 
play an important role and therefore their results are combined. There is, however, a sig- 
nificant overlap between the events selected by the H'^H"— >qqqq and H+H^^AW+*AW~* 
selections, and the events selected by the H+H~— ^qqri/',- and H+H^— ^AW='=*r^',- selections. 
Therefore, an automatic procedure is implemented to switch off the less sensitive of the over- 
lapping channels, based on the calculation of the expected limit assuming no signal. In general 
the fermionic channels are used close to the [m^±,mA] diagonal and for low tan/3, and the 
searches for H='=— >AW^* are crucial for low values of niA and high values of tan /3. 

The confidence 1 — CL^ is calculated combining the 6 j + ^ and 8j searches assuming SM 
branching ratios [42] for the W^* decay and is shown on Figure 9(a). The largest deviation 
1 — CLft = 0.0009 corresponding to 3.1a is reached at mH±=45 GeV and mA=44.9 GeV 
(in the middle of the very narrow lightest strip from mH±=40 GeV to 50 GeV at the 
{mii±,mA) diagonal). However, the mean background shift on the [mjj±,mA] plane amounts 
only to 1.1(7. 1 — CLh for the 4j -|- r channel is shown in Figure 9(b). The largest devi- 
ation 1 — CLj, = 0.013 corresponding to 2.2a" appears for low charged Higgs-boson masses 
(mH±=40 GeV, mA=21 GeV), reflecting the excess of events in the v^=189 GeV search. The 
mean background shift for this channel is O.Scr. Note that the results shown in Figures 9(a-b) 
are model-independent. 

When all channels are combined within 2HDM(I), a few hot spots with a significance 
above 2a survive. This is illustrated in Figures 9(c-d), where the combined results are plotted 
for tan /3=10 and 100. For tan (3=10, the largest excess 1 — CLb = 0.01 corresponding to 2.3a 
is found at mjj±=55 GeV and mA=35 GeV (just before switching from the bosonic to the 
fermionic channels). 

As mentioned previously, the H='=— s-AW='=* decay becomes dominant if the A boson is 
sufficiently light. The smaller tan /3 is, the smaller mA should be. This is clearly seen from the 
structure of the result on Figures 9(c-d): for tan/3=10, the bosonic decay becomes dominant 
at ruA ~ mH± — 18 GeV, while for tan/3=100, it dominates already at uia ~ mH±— 6 GeV. 

The model-independent limits on the charged Higgs-boson production cross section rela- 
tive to the 2HDM prediction are presented in Figure 10(a) for the H'*"H~^AW+*AW~* and in 
Figure 10(b) for H^H^^AW^*r^'T- searches, with the only assumption that W^* decays with 
SM branching ratios. The results combining all channels using 2HDM(I) branching ratios are 
shown on Figure 11 for different choices of tan/3. 



22 



OPAL: H"^ H" ^ AW* AW* OPAL: H"^ H"^AW*tv 




inH± (GeV) nijji (GeV) 



OPAL: H", 2HDM type I, tan|3=10 OPAL: H"^ H", 2HDM type I, tanp=100 




inH± (GeV) nijji (GeV) 



Figure 9: The conEdence 1 - CLh on the [mH±,mA] plane (a) for H+H~^AW+*AW"* 
combining the results of the 6j + £ and 8j searches, and (h) for H"'"H~^AW^*t^'t from the 
4j + T search. The combined results in 2HDM(I) for (c) tan P=10 and (d) tan P=100 are also 
shown. The significance values corresponding to the different shadings are shown by the bars 
at the right. 
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Figure 10: The 95% CL upper limits on the production cross section times relevant and A 
boson decay branching ratios relative to the 2HDM prediction on the [?n.jj± , itia] plane for the 
process (a) H+H-^AW+*AW-* and (b) H+H-^AW=^*rz^^. The plotted curves are isolines 
along which the limit is equal to the number indicated. The expected limits are given for (a) 
BR(H±^AW±*)^ • BR(^^bb)^ = 1 and (b) BR(H±^tz^^) • BR(H±^AW±*) • BR(yl^bb) = 
0.25 corresponding to the maximal value in 2HDM(I). 



For intermediate and large values of tan /3, an unexcluded region appears parallel to the 
[mii±,mA] diagonal both in the observed and expected limits in Figure 11 (see the full and 
dotted black curves corresponding to a cross-section ratio of 1.0). This can be understood 
by making two observations: Close to the diagonal, the eight- or six-jet structure of the 
H^^AW^* signal events becomes less pronounced and the final state turns four-jet-like. As 
the selection variables for the signal and the background become similar, the likelihood cut 
removes more signal events resulting in a drop in efficiency and therefore in sensitivity. This 
decreased sensitivity at the [mH±,"iA] diagonal is clearly visible in Figure 10(a). On the 
other hand, the rate of the fermionic events decreases in 2HDM(I) by the distance from the 
diagonal and by tan P, so the fermionic searches also lose their sensitivity. 

To further study the behavior of the unexcluded regions, 90%, 95% and 99% CL excluded 
areas are shown in Figure 12 for different choices of tan/3. At the 90% CL, the small unex- 
cluded islands at mA=12— 15 GeV disappear and the observed limit (as expected) is set at 
the transition region discussed in the previous paragraph. 

Due to the excess of events in the H"'"H~^AW"''*AW~* searches in the year 1999 data, 
the observed limit is lower than the expectation in all regions where the H^— >AW^* decay 
dominates. 

Finally our results are presented independent of tan/3 in Figure 13, and the limits on 
the charged Higgs-boson mass are summarized in Table 9. The absolute lower limit on the 
charged Higgs-boson mass is 56.5 GeV for tan/? < 100 and 12 GeV < mA < n-Ri) which 
should be compared to an expectation of 71.0 GeV. The observed limit is set by tan/?=3.5. 
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OPAL: H"^ H", 2HDM type I, tan (3=0.1 
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Figure 11: The 95% CL upper limits on the production cross section in 2HDM(I) relative to 
the theoretical prediction on the [mii±,mA] plane for different choices of tan/?; (a) 0.1, (b) 
1.0, (c) 10.0 and (d) 100.0. The plotted curves are isolines along which the limit is equal to 
the number indicated. 
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Figure 12: Excluded areas at 90%, 95% and 99% CL on the [m^±,mA] plane in 2HDM(I) for 
different clioices of tan 0: (a) 0.1, (b) 1.0, (c) 10.0 and (d) 100.0. 
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mH±(GeV) 

Figure 13: Excluded areas in 2HDM(I) on the [mH±,mA] plane independent of tan f3 at 95% 
CL. The weakest overall mass limit is defined by the tan (3=3.5 exclusion, which is also shown. 

where a small unexcluded island is present around mH±~60 GeV and mA~12 GeV (also 
shown on Figure 13). This unexcluded island is no longer present at 90% CL where the mass 
limit improves to 66.8 GeV. 

For toa>15 GeV, the tan /3-independent lower limit on the charged Higgs-boson mass at 
95% CL is 64.8 GeV with 71.0 GeV expected. Both the observed and expected limits are 
found in the transition region where the bosonic and fermionic channels have comparable 
sensitivities. The 6 GeV difference is due to the excess observed in the H"'"H~^AW''~*AW~* 
search. On the (mH±, toa) diagonal, the H^— >AW^* decay becomes kinematically suppressed, 
and the 2HDM(II) result for BR(H±^Ti/r)=0.65 of 81.6 GeV (see Figure 8 and Table 8) is 
reproduced. 

7 Summary 

A search is performed for the pair production of charged Higgs bosons in electron-positron 
collisions at LEP2, considering the decays H'^^tUt, qq and AW='=*. No signal is observed. 
The results are interpreted in the framework of Two-Higgs-Doublct Models. 

In 2HDM(II), required by the minimal super symmetric extension of the SM, charged Higgs 
bosons are excluded up to a mass of 76.6 GeV (with an expected hmit of 76.8 GeV) when 
BRiU^^TUr) + BR(H=^^qq) = 1 is assumed. See Figure 8 and Table 8 for BR(H=^^rz/^)- 
dependent limits. 

In the more exotic 2HDM(I), where fermionic decays can be suppressed and H^^AW^* 
can become dominant, a tan /3-independent lower mass limit of 56.5 GeV is observed for 
mA> 12 GeV (with an expected limit of 71.0 GeV) due to a non-signal-like excess observed 
at ^=192-202 GeV in the H+H"^AW+*AW~* search, discussed in Section 4. For mA> 
15 GeV, the observed limit improves to mH±> 64.8 GeV (with an expected limit of 71.0 
GeV). Figure 13 shows the excluded areas in the [mH±, wa] plane and Table 9 reports selected 
numerical results. 
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tan/3 


ruA 


limit on mH± (GeV) 






observed 


expected 


< 100 


12 GeV< m-A < ?n-H± 


56.5 (3.5) 


71.0 (1.0) 




mA=12 GeV 


56.5 (3.5) 


71.0 (1.0) 




mA=TOH± /2 


66.0 (3.5) 


73.5 (1.5) 




"lA > JTiHi — 10 GeV 


64.8 (100) 


71.7 (100) 




mA > JTijji — 5 GeV 


80.5 (100) 


78.0 (100) 


0.1 


12 GeV< m-A < "t-h* 


81.6 


81.1 




mA=12 GeV 


81.6 


81.1 




mp=mY{±/2 


81.8 


81.5 




ruA > ?TiH± — 10 GeV 


81.9 


81.6 




TTiA > Tnii± — 5 GeV 


81.9 


81.6 


1 


12 GeV< niA < ^h± 


66.4 


71.0 




mA=12 GeV 


66.4 


71.0 




mA=mii± /2 


78.3 


76.6 




mA > 'm}i± — 10 GeV 


81.9 


81.5 




mA > ?Wh± ~ 5 GeV 


81.9 


81.6 


10 


12 GeV< niA < 'tih* 


65.8 


73.9 




mA=12 GeV 


69.0 


82.5 




mA=?TiH± /2 


86.7 


89.8 




mA > JTT-Hi ~ 10 GeV 


81.5 


80.4 




niA > TiRi ~ 5 GeV 


81.9 


81.5 


100 


12 GeV< mA < 


64.8 


71.7 




mA=12 GeV 


69.3 


82.7 




mA=mH± /2 


87.1 


89.9 




niA > — 10 GeV 


64.8 


71.7 




TTiA > TiRi ~ 5 GeV 


80.5 


78.0 



Table 9: Lower mass limits for the charged Higgs boson in 2HDM(I). For the tan f3< 100 
results, the tan /? value at which the limit is set is indicated in parenthesis. For any tan /3 
value, an extrapolation of the exclusion hmits to mH± =mA gives the result quoted in Table 8 
for BR{Il^^TUr)=0.65. 
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